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Abstract 


According to the theory of quantum electrodynamics, photon-photon scattering can take 
place via exchange of virtual electron-positron pairs. Effectively, the interaction can be 
formulated in terms of non-linear corrections to Maxwell’s equations, and hence may be 
treated by classical non-linear electrodynamics. Due to the strong electromagnetic fields 
needed to reach any noticeable effect, photon-photon scattering has not yet been observed 
experimentally, but recent improvements in laser technology have increased the possibility 
of direct detection. A verification of the phenomena would be of great scientific value as 
a confirmation of quantum electrodynamics. 

In this thesis the possibility of direct detection of elastic photon-photon scattering 
through four-wave coupling is investigated, both for current and future systems. It is 
shown how three colliding laser pulses satisfying certain matching conditions, can generate 
scattered radiation in a fourth resonant direction. The interaction is modeled, and the 
number of scattered photons is estimated for optimized configurations of incoming pulses, 
both for the Vulcan laser at the Rutherford Appleton Laboratory as well as for the planned 
XFEL at DESY. An experiment using the Vulcan laser is predicted to produce 22 scattered 
photons per shot, an amount which is definitely detectable. The usefulness of the XFEL 
highly depends on future developments in laser technology, but a realistic estimation gives 
a result of 30 000 scattered photons per second. 
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0.1 Preface 


This master thesis originates from a suggestion set out by my supervisor Mattias Marklund 
in springtime 2004. Although some preparations started in August, the main work begun 
in late October the same year, and has been carried out within the Non-linear physics 
group at the Department of Physics, Umea University. 

I became interested in QED-effects thanks to my participation in a graduate course 
on QFT, which provided me with a basic understanding of photon-photon scattering. 
However, the calculations in this work are mainly performed within the framework of 
classical non-linear electrodynamics, and that is why I initially had to put some efforts in 
learning about non-linear interactions. Thereafter a time of experimentation with different 
configurations and models followed, which finally lead to the results presented here. 

Chapter 1 is an introduction containing a general description of photon-photon scatter¬ 
ing and the purpose and method of this work. There is also a survey of the unit systems 
and conventions used later on. 

In Chapter 2 a short overview of QED is presented, and the lowest order Eeynman 
diagram corresponding to photon-photon scattering is pictured. 

The effective field theory outlined in Chapter 3 allows us to express the QED-effects 
giving rise to photon-photon interactions in terms of non-linear corrections to Maxwell’s 
equations of classical electrodynamics. 

Chapter 4 covers the theory of non-linear interactions needed for the calculations in 
this work, including coupling equations and matching conditions. 

In Chapter 5 non-linear terms are calculated for two dimensional as well as three 
dimensional configurations. 

In Chapter 6 the scattered intensity is determined for and compared between different 
models of the interaction. 

Numerical values are inserted in Chapter 7. The number of scattered photons for 
optimal configurations is estimated and discussed, both for the existing Vulcan laser at 
the Rutherford Appleton Laboratory and for the planned XFEL at DESY. 
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Chapter 1 

Introduction 


In classical theory, electromagnetic waves propagating in vacuum cannot, due to the linear¬ 
ity of Maxwell’s vacuum equations, interact with each other, and colliding light will not 
give rise to any scattering. However, the theory of quantum electrodynamics (QED) opens 
for the possibility of photon-photon scattering via virtual electron-positron pairs in the 
vacuum. 

It is possible to formulate an effective field theory containing only the electromagnetic 
fields, valid for field strengths below the critical field 10^® V/m and wavelengths larger 
than the Compton wavelength 10“^^ m, in which the interaction shows up as non-linear 
corrections to Maxwell’s equations. Hence this pure QED-effect may instead be treated 
with classical non-linear electrodynamics. 

In order to reach any noticeable effects, the fields need to be quite strong, and this is 
the main reason why photon-photon scattering has not yet been detected experimentally. 
Recent improvements in laser technology however make the situation look more promising. 

1.1 Purpose and methods 

1.1.1 Purpose 

The aim of this thesis is to investigate, using parameters for current and future systems, 
the possibility of direct detection of elastic photon-photon scattering. Since it seems like 
today’s technology is just on the edge to make it work, an important task is to carefully 
optimize the experiments. 

Photon-photon scattering is a phenomena signihcant for QED, and direct detection 
would be of great scientific value as a confirmation of the theory, comparable to the 
detection of the Casimir effect [2]. Hence the purpose of this thesis is in first place of 
fundamental nature. Of course future fields of application may arise from the knowledge 
of how to experimentally produce and control this scattering. 

1.1.2 Earlier detection schemes 

To simply collide laser beams has not really been a realistic way to detect photon-photon 
scattering before the recent development of petawatt lasers, but during the last decades 
a number of other detection schemes, based on imposing boundary conditions, have been 
proposed. 
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The correction terms vanish in the limit of parallel propagating free waves, but if 
instead waves trapped between regions of conducting material are used, there will be a 
non-zero effect. One suggestion is detection by wave mixing inside a laser produced vacuum 
channel in an overdense plasma m- Another uses the idea of non-linear excitations of 
new modes inside waveguides p]. However, there are some problems with this approach 
related to the properties of the conducting material, and it may be hard to separate the 
disturbances it gives rise to from the photon-photon scattering effects. In order to rule 
out any other scattering processes it is also important that the vacuum is almost perfect, 
which can be technically hard to achieve here. 

The ultimate detection scheme would be when there is no disturbing matter at all close 
to the interaction region, and that is the case investigated in this thesis. 

1.1.3 Method of this work 

The method of this work is based on the simple idea of colliding laser pulses in vacuum, 
and trying to get some scattered photons out of it. To increase our chances we will send in 
three carefully chosen laser pulses satisfying certain matching conditions, such that some 
scattering is, due to the non-linearities, generated in a fourth direction, i.e. we will deal 
with four-wave coupling. 

This scheme obviously does not suffer from any problems with properties of guiding 
material. It also mirrors the common picture of scattering as bouncing particles in a nice 
way: Instead of just exciting a new mode inside some guiding structure, we can really see 
that the light changes direction. 

In order to get as nice result as possible, we will try to find an optimal configuration 
of incoming beams, and the resulting number of scattered photons will be calculated for 
two interesting lasers: The Vulcan laser and the XFEL at DESY. 

1.2 Units and conventions 

1.2.1 Unit systems 

This work deals with three different unit systems. Instead of sticking to a single one, we 
will use units appropriate for the situation we are considering. 

QED is with advantage formulated in natural units, where the vacuum speed of light 
and Planck’s constant simply are put as c = h =1. Hence this is the unit system chosen 
for Chapter 2. 

Maxwell’s equations are especially nice to work with in Gaussian units, where the 
electric and magnetic fields are of the same dimensions. Gaussian units are used from 
Ghapter 3 up to and including Eq. (6.12), after which we go over to Sl-units to make the 
insertion of numerical values somewhat simpler. More information about how Gaussian 
and Sl-units are related, and how they easily may be translated into each other, is for 
example to be found in Jackson [S]. 
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1.2.2 Conventions 

Some important conventions and notations used in this thesis are lined up below; 

• Four-vector components are labeled by Greek indices, running from 0 to 3. 

• The metric signature is (-1-), i.e. the non-vanishing components of the flat 

spacetime metric tensor are given by r/oo = l, ??ii= ?? 22 = ?/33 = —1- 

• The completely antisymmetric symbol is dehned to equal -|-1 for (/r, i/, cr, p) 

an even permutation of (0,1, 2, 3), —1 for an odd permutation, and 0 otherwise. 

• When writing scalar products of four-vectors the indices are sometimes suppressed, 
e.g. px = pPx^. The spacetime four-vector x^ is usually only written as x. 

• The four-gradient may be written like or simply just as 

• The commutator of two operators A and B is denoted by [A, B] = AB — BA, while 
the anticommutator is written as [A, B]_^_ = AB -|- BA. 

• Complex conjugation is denoted with a star (*) and Hermitian conjugation with a 
dagger (^^). Vectors with complex components are equipped with a tilde (~). 
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Chapter 2 

QED: A short overview 


This chapter is intended to give a short introduction to quantum field theory (QFT) in 
general, and quantum electrodynamics (QED) in specific. Although the main calculations 
in this thesis are carried out within the framework of classical non-linear electrodynamics, 
a brief introduction of the underlying theory will be outlined here. After all, it is a process 
predicted by QED which is considered, and for the reader not familiar with this theory a 
quick overview should be in its place. However, anyone who wishes to skip this part and 
directly proceed to Chapter 3, may do so without getting into problems. 

It must be pointed out that the treatment given below is anything but complete. 
Innumerable things have been left out, both details as well as fundamentals. More pro¬ 
found texts are for instance given by Mandl&Shaw m and Zee m Note that natural 
units {c = h =1) are used in this chapter. 

2.1 Quantum field theory 

2.1.1 Foundation and motivation 

At the turn of the 20th century some remarkable progresses were made in the science 
of physics. In 1900, Planck’s explanation of the black-body radiation provided the basic 
ideas that would finally lead to the formulation of quantum theory. As an attempt to fit 
theory to experimental data, he proposed that the process of emission and absorption of 
radiation by atoms occurs discontinuously in quanta. Einstein went even a step longer 
when he in order to explain the photoelectric effect, concluded that the electromagnetic 
field itself consists of quanta, so called photons. About the same time, in 1905, Einstein 
also formulated his famous special theory of relativity (SR), forever changing our notions 
of space and time. 

Both non-relativistic quantum mechanics (i.e.the quantum theory of classical mechan¬ 
ics) and special relativity have turned out to be very accurate within their respective 
domain of validity. But how to describe a particle moving at high speed, where both 
quantum theory and special relativity seem to be needed? The Schrodinger equation deter¬ 
mining the time evolution in non-relativistic quantum mechanics is not Lorentz covariant, 
and hence needs to be modified in order to be consistent with SR. Several candidates of 
covariant relativistic wave equations describing different kinds of particles exist, but suffer 
from difficulties in their single particle interpretation. 

By quantizing the fields obeying those equations, and thereby going over from single 


particles to a many particle theory, one can get around these problems in a smooth way. 
Anyway, when dealing with the interaction between radiation and atoms in quantum 
mechanics, the Fourier components of the electromagnetic field are often subjected to a 
quantization, giving rise to the photons. By symmetry reasons it then seems natural to 
quantize also all other particles. The many particle view allows for the creation and anni¬ 
hilation of different particles, something unachievable in the single particle interpretation. 

We have here given a very rough overview of the utility of a quantum field theory, but 
the adequateness of it must of course utterly be determined by experiments. Let us now 
have a look at how the quantization of fields is performed. 

2.1.2 Formulations 

There exist two main approaches to quantum field theory: The canonical formulation and 
the path integral formulation. These formalisms often appear to be complementary, and 
a problem hard to solve using one of them may be a lot simpler in the other, and vice 
versa. A characteristic difference between the approaches is the presence of operators in 
one but not the other. 

We will for both formalisms outline the quantization procedure of any classical field 
theory derivable from an action integral 

S{T.) = , ( 2 . 1 ) 

by demanding the variation SS(T,) = 0 for any arbitrary spacetime region S when variating 
the fields (j)r, r =1,..., N, needed to specify the system. That is, the equations of motion 
of the fields can be derived from the Lagrangian density C{(f)r,(f)r,fj.) by using the Euler- 
Lagrange equations 

d f dC 

In the canonical formulation of quantum field theory the usual canonical quantiza¬ 
tion procedure is generalized to work on fields as well. In Heisenberg’s approach to 
non-relativistic quantum mechanics, a system with configuration space variables qj{t) is 
quantized by identifying the canonical momenta Pj{t), and then turning those conjugate 
coordinates and momenta into Heisenberg operators subjected to the commutation rela¬ 
tions 

[qj{t) ,Pk{t)] = iSjk 
[Qjit) ,qkit)] = \pjit) ,Pkit)] = 0 

When dealing with fields, the dynamical quantities corresponding to the qj{t) are just 
the field amplitudes (/>,.(x, t) at each point in space. Hence the system has a continuously 
infinite number of degrees of freedom. By going over from the discrete to the continuous 
system, and by defining the fields conjugate as 7r,.(x, t) = dC/dipr, the corresponding 
commutation relations which have to be imposed on the operator versions of the fields 
and their conjugate in order to quantize the theory, are seen to be 

[(?!),.(x, t ), 7rs(x', t)] = iSrsS{x - x') 

[(/)^(x, t) , ()>s(x', t)] = [7r^(x, t) , TTsi-X.', t)] = 0 




dcpr 


= 0 . 


( 2 . 2 ) 
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The path integral formulation of quantum field theory is an extension of Feynman’s 
approach to quantum mechanics; The amplitude Z of propagation from a point qi to 
a point qp is equal to the sum of for all possible paths q between the points in 

configuration space. Here, 5(g) stands for the classical action of the path q. By denoting 
the summation over all paths starting at qi and ending at qp by an integral over Dg, we 
may write the amplitude as 


Z = j Dq . (2.5) 

Of course the summation over all paths is not a trivial thing to get a grasp of, but we will 
not go deeper into any mathematical details in this short review. 

The extension to fields is now quite obvious: The amplitude Z for a system initially 
in state I to finally end up in the state F is given by 

Z = J = J , ( 2 . 6 ) 

where Dcj) tells us to sum over all possible field configurations smoothly connecting the 
initial and final states, and 5((/>) = / d‘^xC{(j)) is the action corresponding to the configu¬ 
ration (j) of such type. Note that this summation is over a mindboggling set of possibilities, 
and the only integral that can be solved analytically is the Gaussian, corresponding to the 
free field case. 

The path integral formulation is indeed the computationally more interesting nowa¬ 
days, much thanks to its suitability to treat nonperturbative effects, but some features 
like the appearance of particle-like properties are easier to see in the canonical formalism. 
The Lorentz invariance is obviously built in from the start in the path integral formalism, 
while it may be harder to see in the canonical since the imposed commutation relations 
single out equal times. 

2.2 Quantum electrodynamics 

2.2.1 Quantization of the Dirac field 

We now turn to quantum electrodynamics, the quantum field theory of the electromagnetic 
and the electron-positron fields and their mutual interaction. 

A relativistic wave equation describing spin 1/2 fermions, including the electron- 
positron, with mass m was found by Dirac in 1928. This so called Dirac equation can 
in a representation-free way be written 

- m^{x) = 0, (2.7) 

where 7 ^ are 4x4 matrices satisfying the anticommutation relations 


[7^71+= 27 ^^ 

(2.8) 

as well as the Hermicity conditions 



(2.9) 
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and where ^l>{x) accordingly is a four component spinor wavefunction ^a{x), a =1,... ,4, 
but we choose to suppress all matrix and spinor indices. 

The adjoint field ^^{x) = then satisfies the adjoint Dirac equation 

- m^(x) = 0 , ( 2 . 10 ) 

and Eqs. (2.7) and (2.10) may be derived from the Lagrangian density 

C = ^(x) — m) ■0(x) (2.11) 

by taking ■0 q.(x) and '4>aix) as the independent fields in Euler-Lagrange equations (2.2). 

Unfortunately, the quantization procedures outlined above really only works for bosons, 
and some adjustments have to be made for fermions in order to respect the antisymmetry 
of the wave function. 

The expansion of a canonically quantized bosonic field in a complete set of solutions 
leads to the identification of annihilation and creation operators satisfying certain commuta¬ 
tion relations. Eor spin 1/2 fermions we do a similar expansion in a complete set of 
solutions to the Dirac equation, but now instead impose anticommutation relations on the 
expansion coefficients. 

A complete set of plane wave solutions to the Dirac equation (2.7) can be found from 
consideration of a cubic enclosure of volume V with periodic boundary conditions. For 
each allowed momentum p, with corresponding energy = \/m? + p^, there are four 
independent solutions: 


Ur{p) 


^-ipx 


u^(p) 


o^pX 


r =1,2 


^/u ’ ^/u’ 

where the constant four component spinors Ur(p) and Ur(p) satisfy 

= 0, {'y^^Pfj, + m)vr{p) = 0 

and are chosen to be normalized according to 

pO 

4(p) «r(p) = ?^^(p) ^'r(p) = — • 

The quantized versions of the Dirac field and its adjoint field can be written as 

l'c^(p) ii^(p) -h 4(p) ^'r-(p) , 


V’(x) = 


p,r 


VpO L 


( 2 . 12 ) 


(2.13) 


(2.14) 


(2.15) 


ip{x) = 


p,r 


m 


VpO L 


drip) VriP)e-^^^ +clip) Ur{p)e^^^ 


(2.16) 


where the summations are over all allowed momenta p and both spin states (corresponding 
to r = 1, 2), Vr = 'dr = 47°) and the expansion coefficients are subjected to 


Cr (p), 4 (pO] , = K(p),4(p0 


— drsdpp' 


(2.17) 
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while all the other combinations anticommute. 

It follows from (2.17) that we can interpret Cr, cj. and dr, dj. as absorption and creation 
operators of two different kinds of fermions, and physical properties of these particles 
may be deduced from the constants of motion derivable from the Lagrangian density 
(2.11). The conclusion is that the c- and d-operators corresponds to electrons and positrons 
respectively. For example, with the vacuum state |0) defined by Cr(p) |0) = dr(p) |0) = 0 
for all p and r =1,2, the state dj(p) |0) consists of one positron with momentum p in 
spin state r, while the state c|(p') dj(p) |0) also includes an electron with momentum p' 
in spin state s as well. 

In the path integral formalism, the quantization of fermions can be performed by intro¬ 
ducing so called Grassmann numbers. However, that is nothing which will be investigated 
here, and we choose to be satisfied with the canonical approach. 

2.2.2 Quantization of the electromagnetic field 

Expressed in terms of the four-vector potential A^(x), Maxwell’s equations for the free 
electromagnetic field become (see Section 3.2 for definitions) 

rd^H^(x) - (d^H'^(x)) = 0. (2.18) 

The quantization can unfortunately not be performed using the usual Lagrangian den¬ 
sity (3.23) due to problems with conjugate fields being identically zero, but another better 
suited one is given by 

£ = -^ (d,A^(x)) (d‘'A^(x)) . (2.19) 

However, the field equations this Lagrangian density gives rise to are only equivalent to 
Eq. (2.18) in a Lorentz gauge, i.e. we have to add the subsidiary condition 

a^H^(x) = 0. (2.20) 

On identifying the fields conjugate to A^(x), carrying out the canonical quantization pro¬ 
cedure, and finally expanding the free electromagnetic field in a complete set of solutions to 
the field equations (2.18) in a Lorentz gauge (i.e. to d'^di^A^(x) = 0), one gets the quantized 
version 

= E “'(>") «■"" + . ( 2 ^ 21 ) 

k,r 

with summation over all allowed wave vectors k and four polarization states r = 0,... , 3 
with corresponding polarization vectors e(f, the energy = |k|, and expansion coefficients 
as operators subjected to the commutation relations 

Q,j.(k) , al (k ) — Crdrs^kk' 
ar.(k), as(k') = a|^(k),a|(k') =0 

where Co = and Ci = C 2 = Ca = 1- We are thus lead to the interpretation of ar(k) 
and aj(k) as absorption and creation operators respectively for both transverse photons 
(r =1, 2) and longitudinal photons (r = 3). The minus sign present for scalar photons (r = 0) 
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brings about some problems, but can be handled by imposing the Lorentz gauge (2.20) 
(or really an equivalent condition). This will lead to the conclusion that only transverse 
photons may be real, and hence the number of independent polarizations is reduced to 
two, something familiar from classical electrodynamics. 


2.2.3 The interaction and Feynman diagrams 


After this compressed presentation of the quantized free electron-positron and electro¬ 
magnetic fields, we turn to their mutual interaction, something which allows for annihila¬ 
tion, creation and scattering processes to take place. 

Before doing so, we will rewrite the Lagrangian densities of the free fields a little 
bit. It turns out that for the Lagrangian densities (2.11) and (2.19), the energy of the 
vacuum state is not finite. However, since only differences in energy are observable, we 
can redefine the vacuum energy such that all other energies are measured with respect to 
it. The removing of the infinities can be done by normal ordering of the operators, that 
is moving all annihilation operators to the right of the creation operators, treating them 
as if they all commute (bosons) or anticommute (fermions). Hence we will redefine the 
Lagrangian densities as normal products (denoted by N). 

The form of the interaction Lagrangian density Cj can be found following the procedure 
of incorporating an electromagnetic field into non-relativistic quantum mechanics. There 
are good reasons for this form to be the correct one, but ultimately that is of course up 
to experiments to determine. The total Lagrangian density in QED can be written as 
C = Co + Cj, with free and interaction parts given by 


/:n(x) = N 


tpix) - m) tpix) - ^ (d^Af,(x)) (d‘'A^(x)) 


(2.23) 


and 


C/(x) = N [e'ip(x)j'"Af,(x) V’(a:)] , (2.24) 

where (—e) is the charge of the electron. In the same way, the total Hamiltonian may be 
divided into a free and an interaction part, H = Hq + Hi. 

Due to the non-linearity of the interaction, the coupled field equations become really 
hard to solve, and we have to rely on perturbation theory, which is justified thanks to the 
weak coupling between the free fields. To treat this situation, we have to go over from 
our Heisenberg picture to the interaction picture, where the operators continue to satisfy 
the same equations of motion and the same commutation relations as the free fields in the 
Heisenberg picture, but where the state vectors |4>(t)) carry a time dependence and evolve 
in time in accordance with 

jL !$(«)) = |4,(*)). (2.25) 

In the interaction picture the interaction Hamiltonian, Hi{t) = / d^x7ii{x), also becomes 
time dependent, and in QED it is given through the interaction Hamiltonian density 

Hi{x) = —Ci{x) = N[—eV!(x) 7 ^A^(x)'0(a:)] . (2.26) 

For a scattering process we are interested in the states a long time before {t = —oo) 
and a long time after {t =oo) the interaction takes place. These states are related through 
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Figure 2.1: The lowest order Feynman diagram for electron-positron pair creation. 


the S'-matrix defined by |$(oo)) = oo)). The amplitude for an initial state |i) to 

end up in a final state |/) is correspondingly (/| S\i). By solving Eq. (2.25) iteratively, 
the S'-matrix is identified as the expansion 


OO 


s=i; 

n=0 


HI 

n! 


d^xid^X 2 ... d^XnTiUi (xi) Hi (xa) ...Hi (x„)} (2.27) 


where the integrations are over all space and time. The time-ordered product T orders 
the operators in chronological order with earlier times to the right. 

The S'-matrix expansion may be written as a sum of normal products (Wick’s theorem), 
so when calculating an amplitude {f\ S\i) we only need to pick out the parts that are able 
to first absorb the initial particles, and then create the final ones. When the initial and final 
states are specified by particles with given momenta and spin (polarization), it is found 
that the scattering amplitude can be pictured as a set of so called Feynman diagrams in 
momentum space. Each diagram consists of external lines corresponding to the initial and 
final particles of the process, with internal lines and vertices in between. The internal 
lines may be seen as virtual particles present only during the interaction, and the number 
of vertices corresponds to the perturbation order. Often it is accurate enough to perform 
the calculations to first non-vanishing order, and hence neglect all diagrams with a larger 
number of vertices. Photons are usually pictured by wavy lines, while electrons (positrons) 
are drawn as straight lines with arrows pointing along (against) the time direction, which 
we here choose to be from the left to the right. 

The correspondence between scattering amplitudes and Feynman diagrams may be 
summarized in a set of rules, the Feynman rules for QED. Following those we can in 
principle calculate the amplitude for any scattering process to any given order. We just 
draw all topologically different Feynmann diagrams of current interest, and then via the 
Feynman rules translate those into mathematical expressions. Hence we need not deal with 
the equivalent, but less practical, method of identifying terms directly from the normal 
products following from Wick’s theorem. It may be pointed out that the same Feynman 
rules are obtained by using the path integral formulation. 

As an example, the lowest order Feynman diagrams for electron-positron pair creation 
is shown in Figure 2.1. The diagram represents two initial real photons, an intermediate 
virtual fermion, and a final electron-positron pair. 

When going to higher orders one encounters problems with infinities from divergent 
integrals, threatening to ruin the complete theory. Luckily these infinities can be handled 
through the methods of regularization and renormalization, where the infinities are, loosely 
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Figure 2.2; The lowest order Feynman diagram for photon-photon scattering. 

spoken, put into the relations between bare and physical particles instead. This procedure 
works to all orders of perturbation, and hence QED is called a renormalizable theory. 

It is a general feature, following from the form of the interaction Hamiltonian density 
(2.26), that every vertex is built up by two fermion lines and one photon line. Thus 
we have to add more vertices in order to find a diagram which describes photon-photon 
scattering, that is a diagram with two photons both initially and finally. The lowest order 
photon-photon scattering diagram is pictured in Figure 2.2. It shows that the interaction 
takes place via virtual electron-positron pairs in the intermediate stage. 
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Chapter 3 

Effective field theory 


It may often be useful to relate the properties of the quantum vacuum to properties of a 
medium in classical theory. This can be done by integrating out all high-energy degrees 
of freedom, which are invisible in the classical low-energy domain, leaving us with an 
effective action from which classical equations of motion, effectively incorporating the 
quantum effects, can be derived. QED is a particularly nice candidate for this treatment, 
since it is easy to identify and separate the high energy degrees of freedom, corresponding 
to the massive electrons, from the low-energy massless photons. 

3.1 Effective Lagrangian density 

3.1.1 Quantum vacuum 

The vacuum of quantum field theory differs quite remarkably from the classical view, and 
is sometimes pictured as a stormy sea of fluctuations, with virtual particles constantly 
popping up just to quickly disappear again. Since virtual particles cannot be measured, 
a fact making any discussion about whether they really exist or not meaningless, this is 
of course just a mental picture. The vacuum becomes interesting when we disturb it with 
an external field, and observe how it responds. 

For external electromagnetic fields stronger than the critical field strength of 10^® V/m, 
QED accounts for a significant creation of real electron-positron pairs out of the vacuum. 
This phenomena, as well as the detailed structure of the quantum vacuum in general, is 
carefully treated by Greiner et al. [HI. Photon-photon scattering arises in not so strong 
fields, where the electron-positron pairs only appear as virtual particles in intermediate 
stages. 

3.1.2 Constant electromagnetic field 

The effective Lagrangian density for a vacuum perturbed with a constant external electro¬ 
magnetic field was presented by Heisenberg and Euler in 1936 [7], and is often referred 
to as the Euler-Heisenberg Lagrangian density. The situation is of course an idealization, 
because of the infinite extension of the field, but is a very useful approximation if one 
deals with electromagnetic fields that vary little over an electron Compton wavelength. 
This means that the Euler-Heisenberg Lagrangian density may without problems be used 
for plane electromagnetic waves with wavelengths larger than 10“^^ m. 
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The formal derivation was given by Schwinger in 1951 m, where he used the so 
called proper-time method, a particularly elegant regularization procedure which explicitly 
deals only with gauge-covariantly objects, and thereby conserves formal gauge invariance 
throughout the calculations. The vacuum energy may be represented in terms of an 
integral with respect to a scalar variable which has many similarities with the proper time, 
and thereof the name of the method. The full derivation is beyond the scope of this thesis, 
but a nice review is given by Dittrich&Gies [^. Technically speaking we integrate out the 
electrons and positrons coupled to the external field, i.e. we take care of all the internal 
fermion lines in Feynman diagrams like Figure 2.2. This finally leads to an expression of 
the effective Lagrangian density C for constant external electromagnetic fields, which for 
relatively weak fields can be expanded as 


C, — To + (fT -|-... 




+ .. • , 


(3.1) 


where ^ = fie^/457rm^c^. Observe that the expression is given in Gaussian units, which 
are used from now on. Tq is, as will be shown in the next section, the classical Lagrangian 
density, while bC represents the non-linear QED-correction. Since ^ is very small, the 
correction is negligible for weak fields. Also note the vanishing of bC for parallel propa¬ 
gating waves. We repeat that Eq. (3.1) is valid for field strengths below 10^® V/m and 
wavelengths larger than ICT^^ m, conditions which are met for all electromagnetic waves 
considered in this work. 

A new quasi-classical theory, describing the quantum effects in a classical language, 
may now be derived with the effective Lagrangian density C as starting-point. Both 
classical electrodynamics and QED are characterized by Lorentz and gauge invariance, 
but while the original classical theory is linear, quantum effects introduces non-linear self¬ 
interaction terms. These may be identified with a polarization and magnetization of the 
vacuum, as shown below. 


3.2 Effective field equations 

3.2.1 Maxwell’s equations 

In the case of no free charges or free currents, but including possible polarization P and 
magnetization M, Maxwell’s equations can be written 


V • E = -IttV • P, 

(3.2) 

1 (9B 

V X E-h = 0, 

c at 

(3.3) 

V • B = 0, 

(3.4) 

1 flE 1 ^P 

= dvr-——h dvrV X M. 
c at c ot 

(3.5) 
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From these we can directly derive the wave equations 


E = - c^V^E + c^V (V • E) + dvrc^V (V • P) 




dt"^ 

a^E 


+ c^V X V X E + dvrc^V (V • P) 


'’>4164 


= dvrc 


(3.6) 




52 B 

'W 

d'^B 

dt"^ 


- c^V^B + c^V (V • B) 


+ c^V X V X B 


/ 15E _ „ 

= c^V X (-— + V X B 

\ c ot 

( 1 5P 

= dvrc^V X --— + V X M 
\c ot 


(3.7) 


In (classical) vacuum, where P = M = 0, Maxwell’s equations reduces to 


V-E = 0, 


V X E + -— = 0, 
c ot 


V • B = 0, 


XX XX 15E 

VxB-^ = 0. 

c dt 


These equations can be written in covariant form as 

= 0 , 


(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 


+ PTfM,U — 0, (3.13) 

where the antisymmetric electromagnetic field tensors and F^i, = rf^^rjupF^P are given 
by 


/ 0 

~Ex 

Ey 

-E, \ 

Ex 

0 

-B, 

By 

Ey 

B, 

0 

-Bx 

V 

— By 

Bx 

0 y 


(3.14) 
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0 

Ex 

Ey 

E, \ 


~Ex 

0 

-B, 

By 


— Ey 

B, 

0 

H 

cq 

1 

V 

-E, 

— By 

Bx 

0 ) 


(3.15) 


It can easily be checked that Eq. (3.12) gives us Eqs. (3.8) and (3.11), while Eq. (3.13) 
corresponds to Eqs. (3.9) and (3.10). 

The existence of scalar and vector potentials, and A, defined through 


B = V X A, 

(3.16) 

15A 

(3.17) 

= — _ v$, 

c ot 


is implied by Eqs. (3.9) and (3.10). By defining the four-vector potential as 

A^ = ($,A), (3.18) 


it is straight forward to verify that the field tensor can be expressed as 


F —A —A 

In terms of the potentials, Eq. (3.13) is identically satisfied. 
Further on we notice the identities 


(3.19) 


(3.20) 


where the so called dual field tensor is defined as 


= -e^’^^PPcrp = 


( ° 

cq 

1 

By 

cq 

1 

\ 

Bx 

0 

E, 

-Ey 


By 

-E, 

0 

Ex 


V B^ 

Ey 

—Ex 

0 

/ 


(3.21) 


(3.22) 


This explicitly shows the Lorentz invariance of both and E • B. 

We will now show that the field equations (3.12) can be derived from the Lagrangian 
density 


Cq = — (e^-B^]= - 

^ 8tt\ J levr ^ 


(3.23) 


by treating the components of as the independent fields in the variational principle, 
that is solving Euler-Lagrange equations 


d f dCp \ _ dCp 

dxf^ j dA^ 


(3.24) 
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We have 


dCo 


Stt 


1 dFa0 

levr ^ IGir dA^^^ 

1 d {Ap f^—Agj ^) ^ 

Stt dAiy^f, 


1 

Stt 




Stt ^ ^ Stt ^ ^ 47r 


giving 


and since 


9x^ \ <9^1/,At / 


J_ 

dvr 




StCo 

5A, 


= 0 


(3.25) 


(3.26) 


(3.27) 


we finally end up with Eq. (3.12). 


3.2.2 Effective polarization and magnetization 

Now let us investigate what the equations will look like when we consider the total effective 
Lagrangian density from Eq. (3.1); 


C = Co + SC 



B^] + 


Stt 


1 £ 
_ _p _i_ 

167r + Stt 


(^E'^ - B'^y + 7 (E ■ Bf 

+ l-Fp.FP'^F^pF^P^ . (3.28) 


We have 


d(F^^F^PF„pF^p) _^d(F^pF<-y 


dA 


U,fl 


dA 


-F.pF^P 


ly,}! 

dC 


3271—^F„pF^P = 16 (e^- B^) F‘'p, 
oAy^p V / 


(3.29) 


and 


8(F,,7^i‘-F^,}^-'‘) „„ mB{F„0F,s) 

dAjj^p dAyp dAy^p 


_ ^aP'yS 


9 {Ap^a - Aa,/3) , 9 {As^^ - A^^s) jp 

9A,,p + dA,,p 

{s’^pSP - S'^Sy F,s + {S'^5P - 6;sy F^0 
+ (e^PP'^- d^P^p) F 
= AeP^^PF^pF^pF^P = 32 (E • B) F^p, 


_ ^al3'y& 


F^pF^P 


FapF^P 




F.pF^P 


(3.30) 
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so Eqs. (3.29) and (3.30) together with (3.28) yields 
9 fd{SC)\ C 


dx>^ I dAiy^fj_ I dvr L 


We also have 


2 + 7 (E • B) 

d{5C) 


dA,. 


= 0 , 


(3.31) 


(3.32) 


so by putting Eqs. (3.26), (3.27), (3.31) and (3.32) into Euler-Lagrange equations (3.24), 
but now with the total effective Lagrangian density C = Cq + 5C instead of just Cq, we 
arrive at the effective field equations 

^ [-2 F''^^ - 7 (E • B) 

It is straight forward to check that u = 0 corresponds to 

V-E = -^V - [2 (e2-E2)e + 7(E-B)B| 
and u = 1,2,3 is equivalent to 




(3.33) 

(3.34) 


V X B - [2 (e ^-E + 7(E • B) B 

c dt ^cdt I \ ^ ■ 


+ CV X [-2(e2-E2)b + 7(E-B)E 


(3.35) 


Comparing these results with Eqs. (3.2) and (3.5) reveals that we can express the correc¬ 
tions to Maxwell’s equations, arising from photon-photon scattering, as an effective polar¬ 
ization and magnetization given by 

^ 2 { /r2 d2^ 

dvr L 


P = — 2(e^-E^)E-F7(E-B)B 


M = 


dvr L 


-2(E2-E2)B-h7(E-B)E 


(3.36) 

(3.37) 


3.2.3 Effective wave equations 

Thus the wave equations (3.6) and (3.7) can be written 


dF 


^ E = -c^ I V(iVl) + -^(Ar234) 


1^- 
c dt ' 


dF 


B = cW X (A^234), 


where (Nl) and (E234) are defined as 

(iVl) = -^V- [2 (e2-E2)E-F7(E-B)B 

1 d 




2(e2-E2) E-F7(E-B)B 


CV X [-2 (e^- E^) B 7(E • B)E 


(3.38) 


(3.39) 


(3.40) 


(3.41) 
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Chapter 4 


Non-linear interactions 


The interest in non-linear effects has been steadily increasing during the last four decades, 
and it has been especially important for the development of fields such as plasma physics 


and modern optics. Here the theory will only be worked out to the extent that it covers 
what is needed for coming chapters. 

Non-linear interactions may be classified into weak and strong non-linearities. If the 
variation of the wave amplitude due to non-linearity is slow compared to the harmonic 
oscillation we have a weak non-linearity, while it is said to be strong otherwise. 

There are two main methods for describing weak non-linear wave-wave interactions, 
namely the coherent-phase description and the random-phase description. In the coherent- 
phase description the coherence time is assumed to be much longer than the time of inter¬ 
action, while the opposite assumption is made in the random-phase description. Therefore 
phase effects play an important role in the first but not in the later case. 

In this work weak interactions between laser beams are considered, and thus we are 
interested in weak coherent non-linear wave-wave interactions. A more complete treatment 
of this area can be found in Weiland&Wilhelmsson HZ!. However, we will just restrict 
ourselves to the case of four plane waves, of which one is much weaker than the other 
three, interacting in a small region of spacetime. By carefully choosing wave vectors to 
satisfy certain matching conditions, the weak wave can be driven by the others, that is its 
amplitude increases due to the non-linear interaction. Let us see how this comes about. 

4.1 Four-wave coupling 

4.1.1 Coupling equations 

As starting point we consider four plane waves, with amplitudes allowed to have a weak 
spacetime dependence due to interactions. 



together with a non-linear wave equation of the form 
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(4.2) 


8 8 8 

1=1 m=l n=m 

where for i = 5,..., 8 we have defined Ei = E*_^, k* = —kj_ 4 , uji = —uji-i, and 
are constants small enough to ensure weak non-linearity, while are unit vectors. As 

will be apparent, this is in fact the form of our wave equation of interest. We also assume 

that the wave vectors and frequencies satisfy the matching conditions 

ki -h k 2 = ks k 4 , (4.3) 

CVi -|- U!2 = ^3 + CJ 4 . (4-4) 

Due to the weak non-linearity, space and time derivatives of the amplitudes are much 
smaller than of the harmonic parts. Taking the dispersion relation loj = ckj into account, 
we thus have to lowest non-vanishing order 



and clearly 



where ^ + ^ ' ^) = convective derivative along the wave in 

question. 

From Eqs. (4.1), (4.2), (4.5) and (4.6) we have 



l=l m=ln=m 


Let us multiply both sides of Eq. (4.7) by g“*(^i'r-‘^ii) thereafter average over many 
wavelengths and period times. Due to rapid oscillations, all terms with non-zero exponents 
will then average to zero and be negligible. Thanks to the matching conditions (4.3) and 
(4.4), there is a term that will survive on the right hand side, namely the one including 
the product E^E^E/i. Thus there is no need to ever consider the other terms, an insight 
which definitely makes the calculations easier. When evaluating the right hand side of a 
wave equation like (4.2) we only have to keep track of the resonant terms, that is only 
determine one of the in total 120 different Similar arguments of course works for 
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the other three waves. So by neglecting all non-resonant terms, we end up with wave 
equations looking like 




dt"^ 




dt"^ 






Ei(r, t) = C 1 E 2 E 3 E 4 U 1 , 

(4.8) 

E2(r,t) = C2EIE^E4^2, 

(4.9) 

E3(r5i) = ^ 3 ^ 1 ^ 2 ^ 4 X 13 , 

(4.10) 

E4(r, t) = C4EiE2E^\i4, 

(4.11) 


where Cj are the coupling coefficients, and Uj are unit vectors. 

Next assume that initially one field is much weaker than the others, say for instance 
£14 El, E 2 , E 3 . From Eqs. (4.8)-(4.11) we then see that, for interactions occupying 
small regions of spacetime, the driving of the strong field amplitudes may be neglected, 
strengthening the point of view to treat those fields as unaffected by the interaction. Only 
the driving of the weak field amplitude will be of importance. 


4.1.2 Generated field 

Now consider the case when the initial weak field in fact does not exist, but the strong 
fields Ei(r, t), E 2 (r, t) and E 3 (r, t) still satisfy the matching conditions (4.3) and (4.4) for 
some k 4 . In other words the wave vectors satisfy the equation 

|ki -I- k2 - kal = fci -I- £2 - ^3, (4-12) 

and we make the definitions k 4 = ki -I- k 2 — ka and k 4 = ki + k 2 — k^- Prom the above 
discussion it is clear we expect a wave to be generated in the k 4 -direction, but how about 
directions differing only slightly from it? To clear this out we write the generated field 
Eg(r, t) as 

E 3 (r, *) = \ (E 3 (r, t) + E*(r, t)) , (4.13) 

where Eg(r, t) is the part corresponding to E 4 (r, t) . Following the procedure 

of just keeping resonant terms, we will end up with an equation like 

Eg(r,t) = C4EiE2E^U4e<'^"'^“‘^^*), (4.14) 

having the solution pQ 

1 f g*(k4-r-w4tfl) 

j EiE 2 El\t^ - - - dV, (4.15) 
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where the retarded time is defined as tji = t — R = |r — r'| and integration is over all 
space. The term EiE 2 E^\tii is kept inside the integral since the fields will only occupy 
a limited region of spacetime. The subscript indicates that the term really should be 
evaluated at tn instead of t. For now we model the amplitudes to be constants inside 
some interaction region, and zero outside. This is to be seen as a source generating the 
fourth wave. 

Assume that this source region is constituted out of the volume V' in space, and infinite 
in time. The solution may then be written 

1 f 

Bg{r, t) = -^C^EiE2Elii^ - - - dV. (4.16) 

In the radiation zone, r ^ r', we have 

R = |r — r^l = + (r')^ — 2r • r' 1 — 2 ^ ^ r — r ■ r' . (4-17) 

In the denominator of the integrand in Eq. (4.16) it is enough to just use i? ~ r, but 
since /c 4 may be large, the full expression (4.17) will be used in the exponent, giving us 
the radiation zone approximation 

1 ^ -uJ4,t+k4r-k4X-V') 

Eg(r,t) ^C4E,E2E*3^a / - dV 

47rc^ JV' r 

=(4.18) 
47rrc^ Jv' 

The fact that k 4 = /c 4 k 4 , used in the last line, follows directly from the imposed matching 
conditions (4.12). Note that the solution takes the form of an outgoing spherical wave, 
multiplied by a direction dependent factor peaked for f = k 4 . 

4.2 Matching conditions 

4.2.1 Interpretation 

In order to take advantage of the resonance phenomena discussed in the last section, we 
need to find a configuration of plane waves with wave vectors satisfying the matching 
conditions (4.3) and (4.4), or equivalently 


ki + k2 = ks + k4, 


( 4 . 19 ) 


ki + k2 = k3 + ki- ( 4 . 20 ) 

These conditions correspond to momentum and energy conservations respectively. (Com¬ 
pare with the momentum fik and energy hto of a photon.) 

Let us first discuss some questions that may arise if one is not familiar with this kind 
of wave interactions. The form of the matching conditions screams for the interpretation 
of wave one and two as incoming photons, while wave three and four should represent 
the scattered ones. But in the experiment we will have three incoming waves generating 
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Figure 4.1; Definition of 9j and fj. 

a fourth one, so how does this then make sense? Well, the problem has already been 
translated from QED into classical electrodynamics via the effective Lagrangian, and we 
should in fact no longer have to interpret the waves as particles. We are free to choose 
whatever matching conditions we like to, and as long as the calculations end up with 
something useful, one should in principle be satisfied. However, there is a simple way of 
thinking which may help. We can see the third wave as an additional external contribution 
to the particles scattered in the ka-direction, with the purpose of stimulating the scattering 
in the k 4 -direction, i.e. making the right hand side of Eq. (4.11) larger. Another question 
is why it is not enough with just three waves, two incoming and a third being generated? 
Since the non-linear terms are cubic in the fields, the wave vectors then have to satisfy 
matching conditions like 2 ki -|- k 2 = ka, 2 ki +^ 2=^3 or 2 ki — k 2 = k 3 , 2 A:i — k 2 = k^. 
For these to hold the three vectors need to point along the same line, i.e. we are confined 
to a one dimensional configuration, making the experimental set up and measuring very 
hard. (This can with advantage be compared to the classical elastic two particle collision.) 
Following the same line of reasoning, we can also rule out the choice of the matching 
conditions ki + k 2 + ka = k 4 , ki + k 2 + k 3 = k^. 

4.2.2 Two dimensional configurations 

Putting focus on finding configurations satisfying (4.19) and (4.20), we write the wave 
vectors as 

kj = kj sin 6j cos (pjX + kj sin 6j sin (pj-y + kj cos 6jZ, (4-21) 

where 6 j and (pj are defined in the usual way, as in Figure 4.1. In total we have twelve 
unknown parameters {kj, 0j, (pj) and four equations, (4.19) and (4.20), so given eight 
parameters the remaining four can in principle be determined. 

For a two dimensional configuration, which we choose to be in the x-y-plane so that 
= ^2 = ^3 = 7r/2 ^ 6*4 = 7 r /2 giving 

kj = kj cos (pjX + kj sin cpjy, (4.22) 
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Figure 4.2: Diagrammatic representations of the solutions to the two dimensional matching 
conditions for equal frequencies, Eqs. (4-26), (4-27): (a) 4>i = (f) 3 , (I )2 = 4>a] (b) 4>i = (j )2 + '^, 

(/>3 = <^4 + TT . 


the system is reduced to eight unknown {kj, fj) and three equations 

ki cos + /c 2 cos 4)2 = ^3 cos (^>3 + k^ cos (4.23) 

ki sin + /c 2 sin 4)2 = ^3 sin 4‘i + ^4 sin 4a (4.24) 

+ ^2 = ^3 + ki, (4-25) 

leaving us with five independent parameters. 

As a hrst approach we consider the case when ki = k 2 = k^ = k, corresponding to using 
three laser beams of the same frequency. According to Eq. (4.25) we will then get ki = k, 
so Eqs. (4.23) and (4.24) reduces to 

cos 4l + cos 42 = cos 03 + cos 04 (4.26) 

sin 01 + sin 02 = sin 03 + sin 04, (4.27) 

which are satished for 0i = 03, 02 = 4a and also when 0i = 02 + tt, 03 = 04 + vr. The 
hrst case where two of the incoming waves must have the same wave vector, may be seen 
mentally as a rhombus with sides as in Eigure 4.2(a). The second case with two pairs of 
head-on colliding beams can be pictured as in Figure 4.2(b). As we should expect, and 
as can be seen from the coupling coefficients derived in the next chapter, the last solution 
is indeed the most interesting. When modeling the interaction we will examine the two 
extreme cases 0i = 7r/2, 02 = — 7 r/ 2 , 03 = tt, 04 = 0 and 0i = 03 = vr, 02 = 04 = 0, 
visualized in Figure 4.3. It will be shown that the later will give rise to a larger generated 
intensity, and therefore is the favoured conhguration. 

4.2.3 Practical adjustments 

There is however some problems with this latest choice of angles (as there is also with all 
other alternatives discussed so far, but let us focus on this case now). For the first, the k 4 - 
and k 2 -direction coincides, making it really hard to separate the small amount of scattered 
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Figure 4.3; Two specific configurations satisfying the matching conditions (4-26) and 
(4.27): (a)(^i=7r/2, (^2 = -7r/2, 4>3='^, 04=O;(b) (/>! = (^> 3 = tt, 02 = (/) 4 =O. 


photons from the intense external beam. Secondly, the generated wave will have its peak 
straight into the sources of the two opposite traveling waves, yielding problems with the 
placing of the detector. What we want to do is thus to find a configuration looking almost 
as the original (to keep as much of the advantages of the model as possible) but where the 
k 4 -direction is slightly changed to make sure none of the above problems arise. 

Say we want to displace k 4 at least an angle S from both k 2 and the sources of the 
opposite directed waves. What is needed to obtain such a configuration is a small shifting 
of frequencies. To get a feeling for how this can be done we draw a rectangle in which 
the diagonals represent ki and k 4 , while two of the sides correspond to k 2 and k 3 , all in 
accordance with Figure 4.4. Guided by this picture we now choose 4ii = tt + 26, (/>2 = 6, 

4>3 = TT + 6, (/>4 = 0, ki = ki = k, k 2 = k‘i = k cos 6, and show that with these choices the 

matching conditions (4.23)-(4.25) are satisfied. Since 

— [kl cos 4>l + /C2 cos 02 — k^ cos 03 — ki cos 04] 
k 

= cos (tt + 26) + cos 6 cos 6 — cos 6 cos (tt + 5) — cos 0 
= — cos 26 + 2 cos^ 5 — 1 = sin^ 6 + cos^ (5 — 1 = 0, (4.28) 


and 


— [kl sin 01 + /c 2 sin 02 — k^ sin 03 — ki sin 04 ] 
k 

= sin (tt + 26) + cos 6 sin 6 — cos 6 sin (vr + <5) — sin 0 
= — sin 26 + 2 cos 6 sin <5 = 0 (4.29) 

kl + k 2 — ks — ki = k + k cos 6 — k cos <5 — A: = 0, (4.30) 



Figure 4.4: Picture of the adjusted wave vectors. 
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Figure 4.5: Specific three dimensional configuration. 

this is clearly the case. Consequently we must be able to shift two frequencies from k to 
k cos 6 in order to follow this recipe. 

4.2.4 Three dimensional configurations 

The restriction to two dimensional configurations is not as limited as it first may sound. 
As soon as two of the waves are heads on, the four wave vectors will always lie in a plane, 
which we of course are free to call the x-y-plane. Although one might expect this to 
be the optimal situation even when considering three dimensions, we have to investigate 
if there could be something favourable with choosing a configuration not confined to a 
plane. It turns out that for a general case the calculations of the coupling coefficients are 
quite hard to perform, and therefore we will just work out one specific three dimensional 
configuration. 

What should be most interesting is the extreme case when the incoming waves form 
straight angles to each other, one along each of the cartesian coordinate axises. This takes 
us as far as possible from a plane. Let us now see if we can find such wave vectors satisfying 
(4.19) and (4.20). To start with we choose = 7r/2, (f>i=0, 02 =ttI2, (/)i=7r/2, 6*3=0, 
i.e.the vectors lie along the x-, y- and z-axis respectively has no meaning in this case). 
To specify the situation completely two more parameters need to be given, so why not 
choose /ci = A :2 = k. Left to be determined are four parameters which from imposing the 
matching conditions easily can be found as ks = k/2, k^ = 3A:/2, 04 = 7r/2 + arctan(l/2\/2) 
and (l)A = 7r/4. We now have a specific three dimensional configuration, pictured in Figure 
4.5, ready to be investigated later on. 
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Chapter 5 

Calculations of non-linear terms 


In this chapter we derive expressions for the non-linear terms in the effective wave equation 
(3.38). It will be assumed the wave vectors of the incoming plane waves (ki, k 2 , ks) have 
been chosen in such way that the matching conditions (4.19) and (4.20) are fulfilled for 
some k 4 . Thus we can use the method of only saving resonant terms as discussed in 
Section 4.1. For every specific wave configuration the calculations hence end up with an 
expression like (4.14). However, since we will consider more general cases than just single 
ones, the results will by convenience be written on the form 

Eg(r,t) = , (5.1) 

where the vector G is a geometric factor depending only on the directions of wave vectors 
and polarization vectors. As before the plane waves will be written on complex form such 
that 

Ej{r,t) = ^ + E*(r, t) . (5.2) 

An attempt to work out the calculations for the most general case of three dimensional 
geometry and arbitrary polarization soon becomes quite messy. Therefore we put focus 
on a general two dimensional geometry, and only work out the calculations for a specific 
example in three dimensions. As discussed in Section 4.2, this is really not such a limited 
restriction since the optimal situations seem to be when we have head-on collisions, in 
which the wave vectors always form a plane. But at least we should check if there could 
be some advantages with a truly three dimensional configuration. 

5.1 Two dimensional geometries 

5.1.1 Arbitrary incoming angles and polarizations 

Wave vectors confined to the x-y-plane can be written 

kj = kj cos (pjX + kj sin 

where (pj is defined in accordance with Figure 4.1. 

With intentions of being general, the calculations will be done for arbitrary linear polar¬ 
izations of the plane waves. Let 7 j be the angle between the z-axis and the polarization 


(5.3) 
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Figure 5.1; Definition ofjj. 


vector of the Fij-field, in such way it is positive in clockwise direction when looking in the 
kj-direction, as shown in Figure 5.1. From inspection of the picture it can be seen that 
the electric fields then are given by 

Ej (r, t) = Ej (r, t) [sin jj sin fijic — sin 'jj cos fijy + cos 7 jz] , (5.4) 

where obviously 

Ej{r,t) = i . (5.5) 

To avoid any possible confusion about the notation it should be pointed out that \Ej{r, t) \ = 
|E*(r, t) I is the amplitude of the oscillating function Ej{r,t), and nothing else. 

The magnetic fields may be expressed in terms of the electric ones using (5.3) and 
(5.4) in the well-known (and from Maxwell’s equations easily derivable) relation for plane 
electromagnetic waves 


Bj (r, t) = kj X Ej (r, t) = Ej (r, t) [cos 7 j sin fijk — cos 7 j cos fijy — sin 7 jz] . (5.6) 

We are now ready to work out the right hand side of Eq. (3.38). Due to the weak non¬ 
linearity the amplitude variations are much slower than the harmonic oscillations, and 
hence the space and time derivatives is taken to operate only on the exponential parts. 
This just brings down some factors in front of the different terms in the expression, so we 
can conclude it really takes the form of Eq. (4.2). Using Eqs. (5.4)-(5.6) and just keeping 
the resonant parts, the non-linear expression may in a tedious way be calculated term by 
term, leading to the final result 

Bg{r,t) = (5.7) 




-c^V^ 
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where the geometric factor is given by 
1 f [7 

G 2 d = + 2 I 4 ">'3 s™ (71 + 72 ) - sin 73 cos (71 + 72 ) 


G2d = +-| -cos 73 sin 


X [sin ((/)3 — ^ 4 ) cos (^4 — (sin (/)3 — sin(/) 4 )] sin^ 

■7 

+ - cos 72 sin (71 + 73 ) - sin 72 cos (71 + 73 ) 

X [sin {(f )2 — 4>i) cos (/>4 — (sin (/)2 — sin(/) 4 )] sin^ 


■7 

+ - cos 7 i sin (72 + 73) - sin 71 cos (72 + 73) 


[sin ((^1 — ^ 4 ) cos 04 — (sin 01 — sin 04 )] sin 
+ ^ I ^ cos 73 sin (71 + 72 ) - sin 73 cos (71 + 72 ) 
[sin (03 — 04 ) sin 04 + (cos 03 — cos 04 )] sin 
■ cos 72 sin (71 + 73 ) - sin 72 cos (71 + 73 ) 


02 — 03 


01 — 02 


+ 1 ^- cos 72 sin (71 + 73 ) - sin 72 cos (71 + 73 )J 
X [sin (02 — 04) sin 04 + (cos 02 — cos 04)] sin^ 


■7 

+ - cos 7 i sin (72 + 73) - sin 71 cos (72 + 73) 

X [sin (01 — 04) sin 04 + (cos 0i — cos 04)] sin^ 


f [7 . 

+ <! - sin 73 sin 


X sin 


L 73 sin (71 + 72 ) + cos 73 cos (71 + 72 ) 
2 ( 4>3 - 04 \ .2 f 4>1-4>2 \ 


■7 

- sin 72 sin (71 + 73 ) + cos 72 cos (71 + 73 ) 


+ 1 ^- sin 72 sin (71 + 73 ) + cos 72 co; 

X sin" ( sin" ( 


■7 

+ - sin 71 sin (72 + 73) + cos 71 cos (72 + 73) 

XX ■ 2 (4>l - ■ 2 (^ 

Sin ( ^ j sin ( ^ ) J 


The complete calculations are to be found in Appendix A 
the expression at least satisfies some reasonable conditions 


be found in Appendix A. Here we will just check that 

,.xx„x,xxxx,xUlxx xxxxxxxlU^xxxxx, 


5.1.2 Checks of the result and a special case 

Since we are expecting a plane wave to be formed in the k 4 -direction, a basic property 
should be the orthogonality between the generated E-field and k 4 , that is the geometric 
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factor should satisfy 


G • k4 = 0. (5.9) 

From substituting Eqs. (5.3) and (5.8) into (5.9), and using the identity 

sin ((/); - (/)fc) cos^ (j)k - (sin ()>; - sin cos (pk 
+ sin {4>i - (j)k) sin^ + (cos - cos (j)k) sin 

= sin {4)1 - 4)k) - (sin cos 4)k - cos 4)i sin 4)^) = 0, (5.10) 


the orthogonality relation is in fact seen to hold. 

A particularly nice form of G 2 d appears when all three incoming waves are polarized 
in the z-direction, i.e. 71 = 72 = 73 = 0. It then reduces to 


G2<i,o = < sin 


. 2 f ^3 - (p4\ . 2 f 4>i- 4>2 


sm 


+ sm 


. 2 f 4>2 - 4>4\ . 2 f4>i- (43 


sm 


+ sin 


2 (4>i - 4>4\ . 2 [(^2- 4>3 


sm 


z. 


(5.11) 


We will for this special case of the geometric factor check a relation corresponding to 
energy conservation. Since the energy density of a plane wave in vacuum is proportional 
to the square of the electric field amplitude, such a relation should be given by 


A 

dt 



= 0 . 


(5.12) 


Comparison of Eq. (4.5) and Eq. (5.1) allows us to draw the conclusion 

= 2i4u)iG2d,zEiE2El, (5.13) 

and with that obviously 

rl P* 

= -2iiu)^G2d,zElE*2E3, (5.14) 

where G 2 d,z as usual stands for the modulus of G 2 d,z- To evaluate (5.12) we also need 
the total derivatives of the three other amplitudes. Looking back on what has been done, 
these are quite trivial to find. We just have to make some modifications of Eq. (5.13). 
As an example, to get an expression for the derivative of E 3 we have to change to 
EiE 2 E^ to E 1 E 2 E 2 , and also interchange 4>3 and (/>4 in the geometric factor G 2 d,z- The 
other derivatives are obtained in a similar way. Due to its nice form, G 2 d,z will however 
not be affected by the interchanging of angles (just perform it to see), and thus will be 
the same in each of the expressions. Hence we have 
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+uj2ElElE^Ei - UJ 2 E 1 E 2 EIEI 
+uj^EiE2ElEl - iOsElE^EsE^ 

+UJAE 1 E 2 EIEI - lv^EIE^EsE^ 

= {iOl + UJ2 — U)^ — UJ4) EIE2E3E4 

— {(jJi + 102 ~ ^3 ~ ^^ 4 ) EiE2E^E2 = 0, (5.15) 

where the last equality follows from the matching conditions (4.20). We conclude that 
everything is as it should be concerning energy conservation. 

The verifications of conditions (5.9) and (5.12), the first for general polarizations and 
the second for a special case, indicates that the G 2 d we have found at least has some 
expected properties. This is of course nothing but a hint of being on the right track, and 
for the full proof of (5.8) one must turn to Appendix A. 

5.2 Three dimensional geometries 

5.2.1 Fixed incoming angles and arbitrary polarizations 

As discussed above, it would be interesting to see what the geometric factor G will look 
like in a truly three dimensional situation. We choose to work with a configuration which 
is, so to say, as far from a plane as possible. Such a set of wave vectors satisfying the 
matching conditions (4.19) and (4.20) was found at the end of Section 4.2, and may be 
written as 

ki = kx 
k2 = ky 
ka = |z 
k 4 = A:x + ky 

Although the wave vectors now are fixed, we can still investigate the configuration for 
arbitrary polarizations. The polarization directions of the two waves lying in the x-y-plane 
can be defined exactly as before through the angles 71 and 72 , but the same definition is 
obviously not possible for the third one as it is pointing in the z-direction. Instead we let 
/Is be the angle between the x-axis and the polarization vector of the Fs-field, such that 
when looking in the ks-direction it is positive in the clockwise direction, see Figure 5.2. 
The incoming electric fields will now be 

El(r, t) = Ei{r, t) [- sin 7 iy + cos 71 z] 

E 2 (r, t) = £' 2 (r,t) [sin 72 X + COS 72 Z] >, (5.17) 

E 3 (r,t) = £' 3 (r,t) [cos/^sx + sin/^sy] 

and thus giving us the magnetic fields 

Bi(r,t) = ki X Ei(r,t) = Ei{r,t) [-cos 7 iy - sin 71 z] 

B2(r,t) = k2 X E2(r,t) = E2(r, t) [COS72X - sin 72 z] >. ( 5 . 18 ) 

B 3 (r,t) = ks X E 3 (r,t) = E 3 {r,t) [-sin/^sx + cos/^sy] ^ 
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Figure 5.2; Definition of 


Following the same method as for the two dimensional case, that is derivatives only 
operate on the exponentials and we throw away all non-resonant terms, we finally arrive 
at the non-linear wave equation 

Bg{r,t) = (5.19) 



with the geometric factor 


Gsrf — — 


1 

2 


1 

2 


UAXJ. / 


7 
+ 4 


-F sin 72 - i cos 72^ cos (71 -h fis 
+ Q sin 71 -F ^ cos 71^ sin (72 -F fis 
sin /33 - ^ cos sin (71 -h 7; 
+ cos72 - I sin72^ sin (71 7 

/ 1 -I \ 


I ■ 

+ ( I cos 71 - -sm 7 i 


a 

- cos fis 


cos (72 4- 


7 

+ 4 


-sin/^s^ cos (71 -^72) 

-F Q cos 72 - ^ sin 72^ cos (71 -h / 

+ (-^ cos 7i - ^ sin 71^ sin (72 
7 r / 1 \ 


'cosfis -h -sm 


fisj sin (71 


X 
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+ f Jsin72 


+ Qsin7i 


+ 


+ ^ cos 72 ^ sin (71 + (3^) 
^ cos 71 ^ cos (72 + / 33 ) 
^ (cos (5^ + sin/ 33 ) cos (71 + 72 ) 

(sin 72 + cos 72) cos (71 + 13 ‘i) 

O 

(cos 71 - sin71) sin (72 + ( 3 ^) 

O 

7 ri 

+ ^ 2 ^3 - cos/33) sin (71+72) 

(sin72 - cos 72) sin (71 + ( 3 ^) 

O 

+ ^ (- cos 7i - sin 71) cos (72 + /33) | 


(5.20) 


The calculations are shown in Appendix B, and here we just verify that condition (5.9) is 
satisfied. 

5.2.2 Check of the result and a specific example 

Using Eqs. (5.16) and (5.20) we have 

- 2 G 3 rf • k 4 = Q sin /33 - cos A + ^ cos (3^ 

1 1 

2 

. 1 1 

+ I - sin 72 - - cos 72 + g cos 72 

1 1 1 


Ld 

- sin /Js + ^ cos /^s + ^ sin f3s^ cos (71 + 72 ) 

1 


--sin 72 


+ i sin 72 + ^ cos 72 ^ cos (71 + ( 33 ) 


fl . 1 1 

+ sin 71 + - cos 71 - - 

111 


- T COS 71 


1 1 
--811171 + -COS 71 

1 


- ^ sin 71^ sin {^2 + Ps) 

7 f 1 

+- (-sin/33 - - cos /33 + cos /33 

sin/33 + ^ sin Ps - ^ cos sin (71 + 72) 

7/1 1 . 1 . 

+ - (--COS 72 - - sin72 + - sin72 

cos 72 + ^ sin 72 - cos 72^ sin (71 + P3) 

7 fl 1 1 

+ 4 - g sin 71 + -811171 

cos71 “ ^ "T “ sin71^ cos (72 + P3) = 0 


(5.21) 
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since in each bracket the terms exactly adds up to zero. Hence we have checked the 
orthogonality condition (5.9) for this configuration. 

As a specific example we may consider 71 =72 = /33 = 0, that is with wave one and two 
polarized in the z-direction and wave three along the x-direction. The geometric factor 
then becomes 

13. 13.. 13. , , 

Gsdo = —X-yH-z. (5.22) 
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Chapter 6 

Interaction models 


As discussed in Section 4.1, the interaction region needs to be specified in order to integrate 
out the generated fields from the derived wave equations. Any too complex model leads 
to great difficulties in the computation, and we will have to sacrifice some accuracy to 
be able to keep as much of the calculations analytical as possible. After all, we are just 
heading for a nice estimation of the scattered amount of photons, but we must of course 
watch out for getting too sloppy. 

It seems reasonable, at least as a basic approach, to model the spatial form of the 
incoming laser pulses as parallel epipeds with quadratic cross sections, inside which the 
field amplitudes are constant. Although the pulses often are more of a Gaussian shape in 
reality, this should work quite well as long as the energies of them are correct. It could 
seem more natural to use a geometry of a cylinder instead of a parallel epiped, but our 
choice has the advantage of giving interaction volumes well adapted for carrying out the 
integrations in cartesian coordinates, something which simplifies the calculations a lot. 

6.1 Cubic interaction region 

6.1.1 Formation 

Let us consider the case when the region of interaction takes the shape of a cube. This 
will for instance happen when (pi = 7r/2, (p 2 = —'Kj'l and i ;/)3 = vr in our two dimensional 
configuration discussed earlier, provided the beams collide at a single spot, which we choose 
to be the origin. When the beams pass through the origin, a cubic region containing all 
three waves will be formed, exist for a while, and ultimately disappear. The side of the 
cube will equal the side of the quadratic cross section of the beams, while the time of 
existence is determined by the pulse length. We will neglect any scattering during the 
construction and end phases, and only model the situation by an instantaneously created 
cube which after a while suddenly disappears again. Note that our three dimensional 
special case also may be treated by this model, due to the straight angles between the 
incoming waves. 
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6.1.2 Generated electric field 


Following the theory outlined in Eqs. (4.14)-(4.18), the radiation zone solution to the wave 
equation (5.1) will be 

E (r,t) = dF', (6.1) 

vrr Jv' 

where V is the cubic region centred at the origin. Let us call the side length h. Consider 
the two dimensional case where we have 


^k4 — r^ = (cos (^>4 — cos (f) sin 0 ) x + (sin — sin (f) sin 0 ) y — (cos 6 ) z, ( 6 . 2 ) 

and thus for cp ^ (f> 4 ^ and 0 7 ^ 7r/2 the integral in (6.1) can be evaluated as 


r ifc4^k4-f^-r' 

Jv 


dV 


^ik 4 {cos(p 4 —cos(l)sm 9 )x'^^/ I ^ ^i/c 4 (sin 04 —sin 0 sin 0 )y^^^/ j ^ ^—ik 4 (cos 6 )z' 


^ik 4 (cos ^ 4 —cos (p sin 0)x' 

2 

^ik 4 (sin ^ 4 —sin (p sin 6)y' 

2 

g— 2 /^ 4(008 6)z^ 

iki (cos 04 — cos 0 sin 9) 

6 

2 

iki (sin 04 — sin 0 sin 9) 

h 

2 

—ik 4 (cos 9) 


2 sin 


I (cos (j )4 — cos (p sin 9) 


2 sin 


(sin (/)4 — sin p sin 0 ) 


ki (cos (/>4 — cos (p sin 6 ) 


/c 4 (sin 04 — sin p sin 9) 


2 sin 


/C 4 I (cos 0 ) 


k^ (cos 9) 

Plugging (6.3) into (6.1) gives us 


(6.3) 


Bg{r,t) = 




Sin 


/C 4 I (cos 04 — cos 0 sin 0 ) sin A; 4 | (sin 04 — sin 0 sin 0 ) 


A; 47 rr (cos 04 — cos 0 sin 9) 
/C 4 I (cos 9) 


(sin 04 — sin 0 sin 9) 


sin 


G2dElE2-E3e*(^4r--a.4t)_ 


(cos 9) 

The real generated electric held is hence given by 


(6.4) 


Eg(r,t) = 




sin 


/C 4 I (cos 04 — cos 0 sin 0 ) sin /C 4 I (sin 04 — sin 0 sin 0 ) 


k^nr (cos 04 — cos 0 sin 9) 


(sin 04 — sin 0 sin 9) 


sin 


/C 4 I (cos 9) 


(cos 9) 


G 2 d\Ei\\E 2 \\E 3 \ cos (A: 4 r - cv^t + 5), 


(6.5) 


where \Ej\ are the amplitudes of the incoming waves, and 6 is the total phase from the 
complex amplitudes Ej. 

In the resonant direction, 0 = 04 and 9 = 7 r/ 2 , that is f = k 4 , the integral simply 
becomes 



( 6 . 6 ) 
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which leads to 


(6.7) 

vrr 

and thereby we can express the real field as 

Bg resir, t) = — —G2d\Ei\\E2\\E3\ cos (k^r - uj^t + 5). ( 6 . 8 ) 

vrr 

Since —> a when x ^ 0 it follows from Eq. (6.5) that 

ff^‘2h3 

lim Eg(r,t) = —^G 2 d|Ei||E 2 | jEsI cos (fc 4 r - a; 4 t + <5) = re 5 (r, t), (6.9) 

</>^</>4 vrr 

as expected. 

Plots of the generated field Eg(r,t) show a characteristic interference pattern with a 
maximum in the resonant direction. Such plots will be considered later on in the next 
chapter. 

To summarize, the generated held may be written like 

Eg(r, t) = Eo{r, 6, cp) cos (k^r - + 5) G 2 d, (6.10) 


where the direction dependent amplitude is given by 


Eo{r,e,4>) = 




Sin 


A; 4 | (cos ())4 — cos(/>sin 0 ) sin ^ 4 ! (sin (/)4 — sin(/)sin 0 ) 


k^nr (cos (j)^ — cos cf) sin 0) 


(sin ())4 — sin (j) sin 9) 


sm 


/C 4 I (cos 6) 


(cos 6) 

with the resonant direction peaked value of 


G2d|Ei||E2||^3|, 


Eo,resir,0, ^) = Eo[r,'^,(t>^ =^^G2d\El\\E2\\E3\. 
^ 1 J TTr 


( 6 . 11 ) 


( 6 . 12 ) 


Until now the equations have been expressed in Gaussian units, but the input of nu¬ 
merical values will get somewhat smoother if we translate our result into Sl-units instead. 
This is done by following the prescriptions [H] 


E —> \/ 47 reo E 
e —> ,} e 

V47reo 


(6.13) 


giving the Sl-expressions 

1 8^ sin 


Eoir,9,(j)) = 


ki I (cos (/)4 — cos cj) sin 6) sin k^ | (sin (/14 — sin cj) sin 9) 


dvreo k^Trr (cos (^4 — cos (j) sin 9) 


(sin 04 — sin 0 sin 9) 


sm 


/C 4 I (cos 9) 


(cos 9) 


■G20E1IIE2IIE3I, 


and 


Eo,resir) = -^^^G2d\El\\E2m\. 
dvreo vrr 


(6.1d) 


(6.15) 


do 


























6.1.3 Generated intensity 

The omission of all non-resonant terms causes our calculated generated electric field to 
have the same polarization, independent of the traveling direction. This leads to the 
unpleasant property of Eg(r, t) and r in general not being mutually orthogonal, as they 
should be for plane waves. However, for f ~ k 4 the orthogonality is almost fulfilled, and 
the plane wave intensity expression 

I = ]^ceoEl (6.16) 


will be safe to use for all directions close to the resonant. The generated intensity is 
accordingly 


Ig{r,9,(j)) = ^ceoE^(r,6',(/>) 


16^2 


\ sm 


/C 4 I (cos 4>a — cos (j) sin 6) 


sm 


k 4 1 (sin 04 — sin 0 sin 6) 
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(cos 04 — cos 0 sin 0 )^ 

2 r/c4|(cos0)l /I 


16^2 
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■ 2 
I sm 


'^2d 


ceolEil^ -ceo|E2| 
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(sin 04 — sin 0 sin Oy 


/C 4 I (cos 04 — COS 0 sin 6) 


sm 


I (sin 04 — sin 0 sin 6) 




(cos 04 — cos 0 sin 6y 
sin^ /c4 4 {cos 9) 

^Glahhh, 


(sin 04 — sin 0 sin 9y 


{cos 9y 


(6.17) 


with the maximum at resonance 

c 2 

Ig,res{r) = ^^4^2^4 ~y:^^2dhhh, (6.18) 

where Ij are the intensities of the incoming waves. 

Following the same steps, a similar scattered intensity pattern is obtained for the three 
dimensional special case, with the peaked value in the k 4 -direction given by 


C 2 

Ig,res{r) = ^^ 4 ^ 2 g 4 (6.19) 

Note that the maximum intensity is determined by the power of the incoming waves 
Pj = b^Ij, together with the frequency (or really for the three dimensional case, oj) and 
the geometric factor (G^|^). However, the shape of the interference pattern depends 
on the size of the cubic interaction region, that is it depends on b. For large regions 
the generated intensity is more focused towards the resonant direction, while the fringes 
becomes more prominent for a small b. 


6.1.4 Two dimensions vs three dimensions 

Let us make a comparison between our two dimensional and three dimensional configura¬ 
tions. The expressions (6.18) and (6.19) reveals that the two cases should give scattered 
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intensities of the same order, as long as and do not differ remarkably. As a quick 
check we compare with Gg^g. By taking (/)i=7r/2, (/)2 = —vr/2, (j) 3 = tt and (/> 4 = 0 

it follows from Eq. (5.11) that in this case 


_ 


1 + ^ + ^) -2.25, 


( 6 . 20 ) 


and from Eq. (5.22) we have 



( 6 . 21 ) 


Thus it seems like there really are no advantages with using the truly three dimensional 
configuration. Of course one should try to optimize the squared geometric factors by 
choosing other polarizations, but since they are built up by trigonometric functions there 
is in fact not much hope of finding something extraordinary. A glance at Eq. (5.8) gives 
us the insight that we should try to choose the 'jj such that the terms including the factor 
7/4 play a greater role. This can be done by taking 71 = 0 and 72 = 73 = vr/2, and thereby 
increase the squared geometric factor to G^^ ss 3.75. Inspection of Eq. (5.20) inspires to 
choose 7 i = 0 and 72 = /ds = 7 r/ 2 , giving G|^ 3.88, so the two cases are still of the 

same order. Guided by this, we choose from this point to abandon all three dimensional 
configurations, and put focus on two dimensions only. If this does not seem well-reasoned, 
it is just because we have not considered the most interesting two dimensional case yet. 
Let us do that right now. 


6.2 Time dependent interaction region 

6.2.1 Formation 

In the model discussed above, the side length of the cubic interaction region is, due to the 
chosen geometry, determined by the width b of the incoming laser pulses. The peaked value 
of the generated intensity is directly proportional to 6 ®, i.e.to the square of the interaction 
volume. However, since it is the power of the lasers which is limited, an increase of the 
beam cross section area must be accompanied by a corresponding reduction of intensity, 
keeping the peaked generated intensity unaffected. We may now ask ourselves if it is 
possible to find another configuration, with the same beam parameters used as previous, 
but with a considerably larger volume of interaction? If so, we expect the scattered 
intensity to be higher, at least as long as nothing dramatic at the same time happens with 
the geometric factor. 

Since the pulse length often exceeds the pulse width with at least an order of magnitude, 
it would be nice if the size of the interaction region somehow could depend on it. This 
is indeed the case when the incoming pulses travel along the same line in space, two 
overlapping waves in one direction and a third in the opposite, eventually colliding head- 
on at some point (the origin). We choose to arrange our coordinate system such that 
this one dimensional configuration is described hy 61 = 62 = Os = 7 ^-/ 2 , (j)i = (f> 3 = tt and 
(j )2 = 0. Moreover we assume that all three pulses are of the same kind, and in specific 
ki = k 2 = ks = k. The matching conditions (4.23)-(4.25) are then satisfied for = k, 
64 = 7r/2 and (/> 4 = 0. 
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Now try to imagine what this region of interaction will be like. Since wave one and 
three are overlapping, these can be pictured as a single parallel epiped with length L and 
cross section area 6^, moving to the left along the x-axis. Consequently a similar parallel 
epiped traveling in the opposite direction then corresponds to wave two. As before, the 
interaction volume is defined as the part of space where all three incoming pulses are 
present, i.e. the region where the two colliding parallel epipeds overlap. When passing 
through each other, the pulses give rise to an interaction volume of the form of a parallel 
epiped with width and hight b in the y- and z-directions, but with a time dependent length 
l{t) in the x-direction. Let t = LIc denote the pulse duration, and the time t = 0 to be 
fixed at the moment when the edges of the pulses first meet at the origin. The length 
l[t) will thus increase linearly from /(0) = 0 to its maximum 1{t/2) = L, and then decrease 
linearly down to 1{t) = 0 again. 

6.2.2 Retarded time dependence 

From Eq. (4.15) it is clear that what we really are interested in is how the length depends 
on the retarded time tji, and not on t. As we are about to measure the scattered intensity 
close to the i;i) 4 -direction out in the radiation zone, we investigate what l{tR) will look like 
when the observer is situated far away on the positive x-axis. At the very first moment 
when the incoming pulses overlap, some initially generated radiation will leave the region 
with the speed of light, going in the (/> 4 -direction. But wave two is also traveling with the 
same speed in the same direction, so all radiation generated at later times at the edge of 
it will reach the observer exactly at the same moment as the portion initially generated. 
Since wave one and three both have the opposite velocity of the second, the front of wave 
two will be situated outside the interaction region as soon as it has crossed x = L/2, 
and hence no more radiation is generated from that edge. The first scattered radiation 
reaching the observer thus comes from a parallel epiped with length L/2 situated between 
x = 0 and x = Lj2. The scattering occurring at the end of pulse two is obviously the last 
portion arriving at the observer, and a similar mental picture tells us it comes from an 
equally sized parallel epiped, now located between x=—L/2 and x = 0. In the same way 
it is easy to convince oneself that for any intermediate portion of radiation received by 
the observer, the size of the region it comes from is always the same, it is just moving 
continuously with half the speed of light in the negative x-direction. Hence we conclude 
that l{tji) = L/2 during some time interval of length r, and zero otherwise. 


6.2.3 Generated intensity 

An expression for the generated intensity during this time interval is thus found in the 
same way as for the cubic interaction volume, provided we replace the limits of integration 
x = ±6/2 with x = ±L/4. The generated intensity for this configuration then becomes 


Ig{r,9,(t)) 


16 ^^ 1 /c4'j (cos(/)4 — cos(/)sin0) sin^ A:4| (sin04 — sin(/>sin0) 
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( 6 . 22 ) 


43 



with the limit in the resonant direction 


^g,res 


(r) = 


kib^ 


ii)' 


dvr^c^eQ 


Gljihh. 


(6.23) 


6.2.4 Cube vs parallel epiped 

From Eqs. (6.18) and (6.23) we find the ratio between the peak intensities of the parallel 
epiped region and the cubic region to be 

^gpara,res{f) 
hcube,res{'f') 

For most laser pulses L^b, and hence the configuration giving rise to the parallel epiped 
interaction region seems to be the favourable. Let us just make a quick check to ensure 
nothing terrible happens with the geometric factor. For simplicity we choose the polariza¬ 
tions 7i = 72 = 73 = 0 (the optimization will be considered in the next chapter). Insertion 
of 4>i = 4>3 = TT, (j )2 = 4’4 = ^ ill Eq- (5-11) leads to the result 

4. (6.25) 

By comparing this to the geometric factor for a cubic interaction region (6.20), we see 
that in fact 02 ^ 02 ^ so there is really nothing to worry about. 

6.2.5 Validity after practical adjustments 

As pointed out in Section 4.2, our choice of configuration brings about some practical 
problems concerning the detection of the scattered photons. This can, as discussed, be 
overcome by a shifting of some frequencies and angles, leaving us with a configuration 
slightly different from that investigated above. However, for the small changes we are 
going to make, it will be justified to continue using the results derived in this chapter. 



(6.24) 
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Chapter 7 

Numerical results 


The Vulcan laser at the Rutherford Appleton Laboratory, capable of delivering petawatt 
shots, is an interesting candidate to use for detection of photon-photon scattering. As is 
the planned X-ray free electron laser at DESY, based on a technique quite different from 
ordinary lasers. In this chapter we investigate and discuss the detection possibilities for 
both of these alternatives. 

7.1 Vulcan laser 

7.1.1 Development of high power lasers 

Since the 1960s when the pulsed laser was invented, the peak power has increased by twelve 
orders of magnitude through a succession of leaps [HI IH E] • After a few years the first 
free running lasers had been improved, first through the process of Q-switching and then 
by mode locking, such that the pulse duration could be cut down from nanoseconds to 
picoseconds for the same amount of energy. Consequently the power was pushed a million 
times from kilowatts to gigawatts. Although some minor amplifications were accomplished 
during the following two decades, the high intensities of the ultrashort pulses prevented 
any larger increase of the peak power. At the associated intensities of gigawatts per square 
centimeter, non-linear optical responses in the gain media imposed an intensity limit for 
the laser. This is due to that the material refractive index shows up a linear dependence 
on the intensity at such high values, and thus causes a deformation of the beam. Even if 
some progress were made in lowering the pulse duration down to femtoseconds, no increase 
in power could at the same time be achieved. 

Thanks to the technique of chirped pulse amplification (CPA) demonstrated in the 
mid 1980s, this problem was mainly solved. In CPA one basically takes an ultrashort 
pulse, stretches it 10^-10^ times and thereby reduces the intensity, amplifies it, and finally 
recompresses it back to its original shape again. By doing the amplification during the 
stretched (chirped) intermediate state of low intensity, one gets around the problem with 
non-linear optical effects. The stretching can be done by passing the pulse through a 
dispersive delay line, and thereby separate the spectral components of the wave packet. For 
ultrashort pulses the bandwidth gets large, and hence it is important that the amplification 
is the same over a large range of wavelengths. The compression, performed by a conjugate 
delay line, must be perfectly matched to the stretching in order to restore the shape of 
the pulse. 
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The CPA lead to the development of terawatt lasers, and in last years several petawatt 
upgrade projects have started. We will put our focus on the Vulcan laser at the Rutherford 
Appleton Laboratory, which fired its first petawatt shot as recent as 5th October 2004. 

7.1.2 Laser capacity 

The Vulcan Laser is capable of delivering 423 J in a single 410 fs pulse at wavelength 
1054 nm, generating a power of 1.03 PW. The beam can be focused to an intensity of 
1.06-10^® Wm“^ on target, corresponding to a beam cross section area of 0.97TCr^*^ m“^. 
The repetition rate takes on the moderate value of one shot per hour. 

In reality the shape of the pulse is of course more complex than a parallel epiped of 
constant intensity, but in order to get something out in a fairly simple way we will stick 
to the model in Chapter 6 . Anyway, we are not aiming for any detailed calculation of the 
scattered number of photons, but only a nice estimation. 


7.1.3 Scattered intensity distribution 

To carry out the experiment, three incoming beams are needed. Thus we have to split 
the petawatt pulse into three parts, each with a third of the original energy. This may for 
instance be done with the help of two beam splitters, but the exact procedure will be left 
for the experimentalists. We hence start out with three beams, each of the same spatial 
shape and wavelength as the original, but with energy, power and intensity cut down by 
a factor one third. 

In Chapter 6 the most promising geometric configuration was found to be described 
by the parameters 01 = 02 = 03 = 04 = 7 r/ 2 , (j)i = (j) 3 = tt, cj)2 = 4>4= 0 , ki = k 2 = k 3 = k^. 
On arranging the incoming beams like this, we expect a generated intensity as in the 
expressions (6.22) and (6.23). The pulse length is now given by L = cr~ 120/rm, while the 
width is from the cross section area found to be 6~9.8 /rm. We also have the wave vector 
ki = 27r/A ~ 6.0-10® m“^ and the intensities of the incoming waves Ii = I 2 = I 3 = I ~ 
0.35-10^® Wm“^ (where A and I are the wavelength and intensity of the original beam). 
We are now ready to substitute these numerical values into Eqs. (6.22) and (6.23), and 
find out what the generated intensity will look like at different spots. To get a picture of 
the angular distribution of scattered radiation, the ratio 
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(where we have substituted (/>4 = 0) is in Figure 7.1 plotted as a function of 0 and (/>, for 
the current values of the parameters k^, b and L. In the limit of fixing 0 = tt/2, the ratio 
as a function of the deviation cj) from the resonant direction becomes 
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and is plotted in Figure 7.2, where the interference pattern might be seen somewhat 
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Figure 7.1: Angular distribution of the scattered intensity, Eq. (7.1), for the Vulean laser. 


clearer. Fixing f = 0 instead, gives a dependence on the deviation = (9 — ' 7 r/ 2 ) from 
resonance as 
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(7.3) 


i.e. of exactly the same form as the ratio (7.2), as expected from the symmetry of our 
chosen geometric conhguration. 


7.1.4 Number of scattered photons 

The total generated power Pg is found from integrating the intensity over a spherical shell 
(with radius r) centred at the origin, that is 
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Figure 7.2: Scattered intensity distribution as a function of the deviation f from resonance, 
Eq. (7.2), for the Vulcan laser. 

The form Eq. (7.4) is written on tells us that about the same power would be generated 
by an intensity profile of constant magnitude Ig^res spanning a square with the side of an 
angle a. 

A numerical solution of Eq. (7.5) gives the value a = 0.096 for our case of interest. The 
central interference peak in Figure 7.1 is definitely the most dominating, and integrating 
over just that region does not have any noticeable effects on a. It is clear that all but a 
negligible part of the power is generated into the central peak directions. From Figures 
7.1 and 7.2 we see that a is around half of the peak width, something which seems very 
reasonable. If we had not done the integration, this is probably what we would have 
guessed just from inspection of the plots. 

Putting Eqs. (6.23) and (7.4) together yields 

Pg = j^^^kih^L^a^Gl,hl2h. (7.6) 

The time of interaction equals the pulse duration t = LIc, so the total radiated energy is 
given by 

Ug = PgT = (7.7) 

Each generated photon carries an energy of Uph=hck 4 , and hence the estimated number 
of scattered photons per shot becomes 
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7.1.5 Optimization of the geometric factor 

All there is left to find is a value for the geometric factor G 2 d given by Eq. (5.8). Due to 
the one dimensional structure of the configuration, the optimization is in this case simple 
to perform. Since all waves travel along the x-axis, we can always arrange our coordinate 
system such that the optimal G 2 d lies along the z-axis. In fact, due to the symmetry of the 
configuration there must be geometric factors of exactly the same magnitudes pointing in 
every direction perpendicular to the x-axis. Anyway, we choose to maximize the magnitude 
of the z-component, which with the current angles imposed takes the form 

7 

{G 2 d)z = I sin 73 sin (71 -h 72 ) cos 73 cos (71 -k 72 ) 

7 

-I-- sin 7 i sin (72 + 73) + cos 71 cos (72 -|- 73) 

3 3 

= - sin 73 sin (71 -h 72 ) + - sin 71 sin (72 -h 73 ) 

cos (71 72 - 73) + cos (72 + 73 - 7 i) • (7-9) 

An upper limit is consequently ( 62^)2 = 3.5, which in fact can be achieved by taking 
7 i = 73 = 7’'/2 and 72 = 0. It is easy to check that the x- and y-components in this 
case vanish, as they should. Hence we conclude that we have found an optimal choice of 
polarization angles. 

7.1.6 Final estimations and practical adjustments 

Finally we plug in all the numerical values into Eq. (7.8), and find the estimated number 
of scattered photons per shot to be 22 , an amount which is definitely detectable. 

It may be pointed out that in the earlier discussed case with a cubic interaction region, 
the width of the central peak becomes almost the same as above, and hence a does not 
change very much. Using Eq. (6.24) we conclude that the number of scattered photons is 
reduced by at least two orders of magnitude compared to the above result. 

For practical reasons, the choice of = ^3 = tt and 02 = </'4 = 0 have to be slightly 
altered. We want the detector to be able to pick up all of the intensity generated into the 
central interference peak, i.e. we have to make sure no other beams are present in that 
region. Following the discussion in Section 4.2, such a modified configuration is given by 
4>i= TT + 2S, 02 = <5, 03 = TT -|- S, 04 =0, fei = ^4 = /c, ^2 = ^3 = ^ COS S, where 5 has to 
be chosen as half the peak width. In this case it means (5 = 0.1. Thus we must shift the 
frequencies of beam two and three by a factor of cos 0 . 1 , giving us the new wavelengths 
^2 = ^3 = A/cos 0.1 ~ 1059 nm and wave vectors ^2 = ^3 = 5.9T0® m“^. As before, we leave 
the technicalities of how to perform this shift to the experimentalists. 

Since such small modifications are done, we will continue to use the same interaction 
model, leading to Eq. (7.8). All parameters are the same as earlier, except for which 
has to be recalculated for the new angles. Putting 0i = tt -|- 0.2, 02 = 0.1, 03 = tt -|- 0.1 
and 04= 0 into Eq. (5.8), with 71 = 73 = tt/2 and 72 = 0, gives G 2 d~ 3 . 47 z. Actually this 
brings the number of scattered photons per shot down a bit, but the rounded off value 
will still he N^ 22 . 
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7.1.7 Suggestion of an experiment 

It is time to summarize the situation, and suggest a single experiment in which photon- 
photon scattering may be detected. The recipe goes as follows: Take the petawatt pulse 
of wavelength 1054 nm produced by the Vulcan laser, and split it into three equal parts. 
Shift the wavelengths of the second and third pulses to 1059 nm. Let the pulses enter the 
x-y-plane at the angles (/>i = vr -|- 0 . 2 , (j )2 = 0.1 and (/>3 = tt + 0 . 1 , and with 71 =73 = vr /2 
and 72 = 0 , corresponding to polarization of pulse one and three in the x-y-plane while 
pulse two is polarized in the z-direction. Also synchronize the pulses to reach the origin 
at the same moment. On the x-axis, place a detector which at least spans the intervals 
- 0.1 < {e-TT/2) < 0.1 and - 0.1 < (j) < 0 . 1 . 

For each shot, we expect several scattered photons to be detected. (More precisely 
our model predicts the average number to be 22 , but we have to keep in mind that a lot 
of approximations have been made. What is important is that the amount of generated 
photons definitely is detectable.) 

7.2 X-ray free electron laser 

7.2.1 Description 

Another interesting candidate to use for detection of photon-photon scattering is the 
XFEL (X-ray Free Electron Laser). Such a laser is planned as a part of the design of the 
superconducting electron-positron collider TESLA (TeV-Energy Superconducting Linear 
Accelerator), at DESY (Deutsches Electronen-Synchrotron), Hamburg [13j . The facility 
is scheduled to take up operation in 2012. TESLA mainly consists of a 33 km long tunnel 
containing two linear accelerators, one for the electrons and one for the positrons, meeting 
at a collision site in the middle. The collider is capable of initially 500 GeV total energy, 
but is extendable to 800 GeV. The first low energy part of the electron accelerator can 
also provide an electron beam for the XFEL. 

The XEEL is based on the principle of SASE (Self-Amplified Spontaneous Emission), 
which can be summarized as follows: The high energy electrons coming from the accel¬ 
erator is guided into an undulator, an arrangement of magnets forcing the electrons to 
follow a winding periodic slalom course. Due to these accelerations the electrons emit 
concentrated X-ray flashes, which will overtake the electrons ahead of them. The passing 
by of the photons causes the electrons to organize themselves into thin layers, and at the 
end of the undulator this structure is fully developed. The electrons in a certain layer will 
now emit radiation in a synchronized way, and hence produce high-intensity ultra-short 
X-ray pulses with laser light properties. This is quite different from how a traditional laser 
works, and one of the main advantages is that the wavelength easily may be altered by 
adjusting the electron acceleration. 

The XFEL attracts scientists from many different disciplines including physics, chem¬ 
istry, material science and structural biology. Here we will investigate how well it can be 
used for detection of photon-photon scattering. 

7.2.2 Laser capacity 

According to the predicted photon beam parameters, the XEEL will be able to produce 
24 GW pulses at wavelength 0.1 nm, with 100 fs duration and a pulse width of 110 /rm. 
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Figure 7.3: Scattered intensity distribution as a function of the deviation f from resonance, 
Eq. (7.2), for the original beams of the XFEL. 

The repetition rate will be as high as 30 000 pulses per second. However, the expected 
coherence time is only 0.3 fs, and hence it should be more accurate to model each pulse 
as a stream of subpulses, each with a length equal to the coherence length of 90 nm. The 
repetition rate may then be expressed as 10^ subpulses per second. 

7.2.3 Estimated number of scattered photons 

Guided by the above properties, and on following the same method as for the Vulcan 
laser, we model the incoming subpulses as parallel epipeds with parameter values given by 
L = 90 nm, 6 = 110 /um, = 6.3-10^'^ m“^ and Ii = l 2 = 1^= 0.66TO^® Wm“^. Figure 7.3, a 
plot of the intensity distribution (7.2), guides us to assign the value a =10“®. Using the 
same geometrical factor as in the last section, Eq. (7.8) finally gives us V~10“^^. Taking 
the repetition rate into account, about 10“^® photons will be generated per second, and 
we may have to wait more than 30 million years to detect the first scattered photon. It 
looks like the high frequency and repetition rate of the XFEL are unable to compensate 
for its low power. 

The configuration used here is however not the optimal. Since we are dealing with 
pulses which are wider than long, we should instead consider the case with a cubic inter¬ 
action region, as in Section 6.1. By performing the substitution L/2 ^ 6 in Eqs. (7.1), 
(7.2) and (7.5) we get the correct formulas for the intensity distribution and c? in this 
case. In fact a plot of the intensity distribution as a function of the deviation from the 
resonant direction will look almost exactly like Eigure 7.3, and we may continue to use 
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Figure 7.4: Scattered intensity distribution as a function of the deviation (f from resonance, 
Eq. (7.2), for the focused beams of the XFEL. 

a = 10“^. Hence we can estimate the scattering number for the optimal configuration 
from Eq. (6.24). Inserting ^ = 3.75, and with the other parameter values as above, 
we see that the scattering number is increased by a factor 1.8-10^ relative to the first 
result. Anyway, the scattering rate is still far to low, making it practically impossible to 
detect photon-photon scattering through this experiment. 

Despite the negative results above, there is still some hope left. It may be possible, 
with future technology, to focus the width of the X-rays down to the theoretical diffraction 
limit of the order of one wavelength, i.e. A M- As a realistic value of the pulse width 
achievable in 2012 we take b = 1.1 nm, a compression of the original value by 10® times. 
Correspondingly the intensity of each beam is increased by a factor 10^*^. From Figure 
7.4, where Eq. (7.2) is plotted for these new values, we estimate a = 0.05. The scattered 
number of photons per subpulse predicted by our model, given by Eq. (7.8), becomes 
N ~ 0.003. Hence we arrive at the incredible result of 30 000 scattered photons per 
second. 

7.3 Discussion 

7.3.1 Improvements 

We cannot expect our simple model to predict any exact scattering numbers, but only 
estimate the order of magnitudes. In order to produce more precise numbers and intensity 
distributions to compare with experimental outcomes, we have to make some improve¬ 
ments. The modeling of the laser pulses as parallel epipeds of constant field strength 
serves our purpose of keeping the calculations analytical and simple well, but a more 
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realistic approach would be to use pulses with a Gaussian-like distributed electric field 
strength, both in the radial direction as well as in the direction of propagation. This of 
course complicates the shape of the interaction region, and we have to rely on numerical 
calculations. No Gaussian models will be considered here, but such calculations would be 
an interesting follow up to this work. Moreover, we have not really taken the effects of 
the bandwidth into detailed account. This may especially be important for the XFEL, 
due to its short coherence time, and a deeper investigation would possibly lead to some 
modifications. 

The improvements may however be overshadowed by practical errors, such as diffi¬ 
culties to tune the laser pulses exactly as we want them, or problems with creating a 
sufficiently good vacuum. There is really no meaning in calculating very precise scat¬ 
tering rates, as long as the experimental uncertainty widely exceeds the accuracy of the 
theory. In the first place we are really just asking ourselves whether some scattering can 
be measured at all, and the model used in this work should definitely be good enough to 
answer such a question. Anyway, what could be compared with experimental results is 
how the scattering number depends on the polarizations of the incoming beams. As long 
as the pulses in reality are symmetric around the axis of propagation, the ratios between 
the number of scattered photons for different sets of polarizations are independent of how 
the pulses are modeled. A sequence of experiments, each with different polarizations but 
otherwise the same, would hence work as a test of the results derived in Ghapter 5. 

7.3.2 Other experiments 

Although this thesis has been pervaded by adjustments to suit the Vulcan laser and the 
XFEL well, there should be no big problems to modify the calculations to fit other situa¬ 
tions as well. For instance it may be practically favourable to use different kinds of lasers, 
with different beam parameters, in the same experiment. We may not be able to follow 
exactly the same procedure of finding resonant configurations, and it may not be possible 
to model the interaction region in exactly the same way as here, but the basic approach 
will be the same, and the modifications needed should be quite easy to figure out. 
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Appendix A 

The two dimensional geometric 
factor 


Here follows the detailed calculation of Eq. (5.7). We calculate the right hand side of 
Eq. (3.38) step by step, keeping only the resonant terms; 
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(E • B) — (El + E2 + E3) • (Bi + B2 + B3) 
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, j , ^ ■ 2 f (1^2 

cos 73 cos 03 cos (71 + 72) Sin ' 


\ 2 f — ^3 

+ cos 72 cos 02 cos (71 + 73) sin ' 


+ cos 


7i cos 01 cos (72 + 73) sin^ a;40. 


E 1 E 2 E 3 


( I ^ ■ 2 f (t'l- 4>2 
Sin 73 cos (71 +72) sin 


. 2 f4>i- (t>3 


+ sin 72 cos (71 + 73) sin^ ^ 

I • ^ I \ ■ 2 f ^2 - ft's 

+ sin 7i cos (72 + 73) sin ' 


0(k4'r-u;40- (^^4) 


E 1 E 2 E 3 


■ i ■ ( , ^ . 2 f 4>i - 4>2 

sm 73 sin 03 sin (71 + 72) sin ' 


\ 2 f — ^3 

+ sin 72 sin 02 sin (71 + 73) sin ' 


\ 2 f ^2 — 03 

+ sin 71 sin 0i sin (72 + 73) sin 


+ 


E 1 E 2 -E 3 


sin 73 cos 03 sin (71 + 72) sin^ 


2 / 01 - 02 


2 / 01-03 


+ sin 72 cos 02 sin (71 + 73) sin^ 


I ■ A ■ / I \ ■ 2 f 4'2 - 4>3 

+ sm 71 cos 01 sm (72 + 73) sm ' 


E 1 E 2 E 3 


cos 73 sin (71 +72) sin^ 


2 / 01-02 


2 / 01 - 03 


+ COS 72 sin (71 + 73) sin^ ^ 

I ■ f I ^ ■ 2 f 4’2- 4>3 

+ cos 7i sm (72 + 73) sm ' 


„i(k4-r-a;4t) 


X 


„i(k4-r-a;4t). 


7 (k 4 -r-LJ 4 t) • 


(A. 5 ) 
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(E B)B = 


V • B- 


■E 1 E 2 -E 3 

4 


+ 


E 1 E 2 EI 


cos 73 sin (/)3 sin (71 + 72) sin^ 

+ cos 72 sin (J)2 sin (71 + 73) sin^ 

+ cos 71 sin (/>! sin (72 + 73) sin^ ^ *^2 ^ </> 3 ^ 

'</>! - 02 ^ 


+ 


E 1 E 2 E 3 


cos 73 cos (^>3 sin (71 + 72) sin^ 

+ cos 72 cos 4)2 sin (71 + 73) sin^ 

+ cos 71 cos 4 >i sin (72 + 73 ) sin^ 

■ I . ^ ■ 2 f ^2 

sin 73 sin (71 + 72) sin ( - - - 


2 

(/>! - </’3 ' 
2 

2 I 422 - 4 > 3 ^ 


i(k4-r-aJ4l) 


X 


gi(k4-r-a;4t). 


+ sin 72 sin (71 + 73 ) sin 
+ sin 71 sin (72 + 73 ) sin^ g*(k 4 -r-a; 4 t) - 


E 


E 1 E 2 E* 

- - — -ik 4 [sin 73 (sin 4>3 cos 4 >a — cos 4>3 sin 4 >a) 

X cos (71 + 72 ) sin^ 

+ sin 72 (sin 4>2 cos 4 >a — cos 4>2 sin 4 > 4 ) 

X cos (71 + 73 ) sin^ 2 
+ sin 71 (sin cf)! cos 4 >a — cos 4 >i sin 4 >a) 

X cos (72 + 73) (^i^)] c-C"* 


E 1 E 2 -E 3 . 
2 


i/C 4 


5 (71 + 72 ) sin ((/)3 - 4 ) 4 ) sin^ 


sin 73 cos 

+ sin 72 cos (71 + 73) sin (02 - 4 >a) sin" 
+ sin 71 cos (72 + 73) 


" 01 -4>2 \ 

7 2 ) 
( 01 - 03 ' 


X sin (01 - 04 )sin 2 ( ^2 ^ 03 ^ g*(k 4 -r -<^40 (^. 7 ) 






= — zt(;4 


E 


(A. 8 ) 
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V X 




- 6 A /4 




+ sin72 cos (71 + 73) sin(^>4 sin^ 
+ sin7i cos (72 + 73) sin ^4 sin^ 


+ 


^ 1 -^ 2 -E 3 


■iA:4 


sin 


73 cos (71 +72) cos (/>4 sin^ 


2 

4>i — 4>2' 


2 f — ^3 

+ sin 72 cos (71 + 73 ) cos ^4 sin ( -^- 


^ 02 — 03 ' 


E 1 E 2 EI . 
2 


iA :4 


+ sin 7i cos (72 + 73) cos 04 sin ^ ^ 

( I N / , A \ ■ 2 I 01 ~ 02 

COS "^3 COS I _ /E , I cm I _ 


,i(k4-r-LJ40^ 


(71 + 72) cos (03 - 04) sin^ 

+ cos 72 cos (71 + 73) cos (02 - 04) sin^ 
+ COS71 cos (72 + 73) 


2 J 

01 - 03 ' 


X COS (01 - 04 ) Sin^ (A. 9 ) 


V X [(E • B) E] = 


■iki 


cos 73 sin (71 + 72) sin (j)4^ sin^ 

72 sin (71 + 73) sin 04 sin^ ^ 0L__^ 

02 — 03 ' 


+ COS' 


+ cos 7i sin (72 + 73) sin 04 sin^ 


+M^ik, 


cos 


73 sin (71 + 72) cos (j) 4 ^ sin^ ^01_^ 

2 0 01 “ 03 ' 


7(k4-r-aj4t)^ 


+ 


E 1 E 2 EI . 


2 / 4^1 — 4^3 \ 

+ cos 72 sin (71 + 73) cos 04 sin ( - - -j 

+ cos 7i sin (72 + 73) cos 04 sin^ ^ 02 ^ 03 ^ ^i{k 4 -r-uj 4 t), 

ill + 72) cos (03 - 04 ) sin^ 

2 0 01 ~ 03 ^ 


■iki sin 73 sin 

+ sin72 sin (71 + 73) cos (02 - 04 ) sin' 
+ sin7i sin (72 + 73) 


X cos (01 — 04 ) sin^ g*(^4'r (A. 10 ) 
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V • [(E • B) B] = 




X sin (71 + 72) sin^ 

+ cos 72 (sin (j)2 cos 4>4 — cos (/>2 sin (/)4) 
X sin (71 + 73) sin^ 


\ ^ / 

+ cos 7i (sin (/>! cos 4>4 — cos (/>i sin (/)4) 

2 02 — </>3 \ 1 iCk^-Y-uJit] _ 

. 2 r" 


X sin (72 + 73) sin^ 

E1E2EI • ^ ^ 

-^-^A;4 [cos 73 sin (71 + 72) 


^01 - 02' 


X sin ((/)3 — (^4) sin^ 

+ cos 72 sin (71 +73) 

X sin (02 — 04) sin^ 

+ COS 7i sin (72 +73) 

X sin (01 - 04) sin^ 1^ 02 - 03 ^ gi(k4-r-a;40 


. 2 . 

^01 - 03 ' 
. 2 . 

^ 02 — 03 ' 




[(E • B) B] = -ia;4 [(E • B) B] 


(A.12) 


(iVl) = ^EiE2^iki I 


■ cos 73 sin (71 + 72) - sin 73 cos (71 + 72) 


2 / 01 - 02 


+ 


+ 


X sin (03 — 04) sin 

7 

- cos 72 sin (71 + 73) - sin 72 cos (71 + 73) 


s. ■ lA. X 1 • 2 01 “ 03 

X sin (02 — (pi) sin ' 


- cos 71 sin (72 + 73) - sin 71 cos (72 + 73) 

X sin (01 - 04) sin^ I g0k4-r-<^4t) (^.13) 
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X • (iV234) 


y • (iV234) 


: —^EiE2E^ik4 |sin73 cos (71 + 72) (siiKp^ — sin(/)4) sin^ 
+ sin 72 cos (71 +73) (sin (j)2 — sin (j)4) sin^ 


^ 4>i — 4>2 \ 

^ 2 J 

^4>i - 4>3' 


+ sin7i cos (72 + 73) (sin 01 - sin04) sin^ ^ *^2 ^ 03 ^ | mt) 

(71 + 72) (sin 03 - sin 04) sin^ 


+^^EiE 2 E^ik 4 I cos 73 sin 

+ cos 72 sin (71 + 73) (sin 02 - sin 04) sin^ ^ 01__^ 


: ^EiE2E^ik4 


+ 


+ 


+ cos 7i sin (72 + 73) (sin 0i - sin 04) sin^ ^ 02 ^ 03 ^ | ^i(ki-Y ui 4 t) 
'7 

- cos 73 sin (71 + 72) - sin 73 cos (71 + 72) 

X (sin 03 — sin 04) sin^ ^ 0i ^ 02 ^ 

■7 

- cos 72 sin (71 + 73) - sin 72 cos (71 + 73) 

X (sin 02 — sin 04) sin^ ^ 0i ^ 03 ^ 

■7 

- cos 71 sin (72 + 73) - sin 71 cos (72 + 73) 

X (sin 01 — sin 04) sin^ / 02 ^ 03 \ | g^(k4•r-a;4^) 


(A.14) 


: ^EiE 2 E^ik 4 I sin 73 cos (71 +72) (cos 03 — cos 04) sin 
+ sin 72 cos 
+ sin 71 cos 


^01 — 02' 


^- iKj - \ I A. ' /^/V-ru - T ~t/ - 2 / 

(71 + 73) (cos 02 - COS 04 ) sin^ ^ 01 _^ 

(72 + 73) (cos 01 - cos 04) sin^ ^ 02 ^ 03 ^ I g*(k4-r-a;4t) 


-^^^1^2^3^^4 |cOS 73 sin (71 + 72) (cos 03 - cos 04 ) sin^ ^ 01 ^ 02 ^ 

1 (71 + 73) (cos 02 - cos 04) sin^ ^0L_^ 


+ cos 72 sin I 


+ cos 7i sin (72 + 73) (cos 0i - cos 04) sin^ ( 0H__^^ | g*(k4-r i^40 


-^EiE 2 E^ik 4 I ^ cos 73 sin (71 + 72) - sin 73 cos (71 + 72) 


+ 


/ J. J, 1 • 2 01 “ 02 

X (cos 03 — cos 04) sin I - - - 

7 

- cos 72 sin (71 + 73) - sin 72 cos (71 + 73) 
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2 f 4 > 1 - 4>3 


+ 


X (cos ((>2 — COS (^4) sin' 

7 

- cos 7i sin (72 + 73 ) - sin 71 cos (72 + 73 ) 

X (cos (pi — cos (pi) sin^ (A.15) 


i ■ (Af234) = 


-2^EiE2Eliki |cos73 cos (71 + 72) sin^ 

g.^2 ( (^2-(pz \ 1 gi(k4-r-a.4t) 


+ 


COS 72 cos (71 + ' 


cos 7i cos (72 + 73) sin^ 


sin (71 
+ sin 72 sin (71 + ' 


\ -I J \ -I J i 

|sin73sin (71 + 72) sin^ 

-73).m^(^^)sin=(^^) 

, ■ ■ ! , I ■ 2 f ~ ■ 2 f 4>2 - (p 3 \ 

+ Sin 71 Sin (72+73) Sin ( - - - I sin ( - - - I 

2^EiE2E^iki I ^ sin 73 sin (71 + 72) + cos 73 cos (71 + 72) 

^ ((>3 - 04 ^ ^.^2 ^ 01-02 ^ 

(71 + 73) + cos 72 cos (71 + 73) 

^ 02 - 04 ^ ^.^2 ^ 01-03 ^ 

■7 

- sin 71 sin (72 + 73) + cos 71 cos (72 + 73) 

*^2 - 03 \ 1 gi(k4'r-aj4t) 


+ 


+ 


X sin 


7 

- sin 72 sin 


X sin 


2 f 04 \ . 2 

X sin - sin 


,i(k4-r-LJ4t) 


(A.16) 


V(A^l) = i/c 4 (iVl) (cos 04 X + sin 04 y) (A-17) 


d - - 

— (Ar234) = -iu;4(iV234) (A.18) 

By finally plugging all this into the right hand side of Eq. (3.38), we end up with the 
generated field of Eq. (5.7). 
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Appendix B 

The three dimensional geometric 
factor 


The details of the calculation of Eq. (5.19) are shown here. On keeping only resonant 
terms, we evaluate the right hand side of Eq. (3.38) step by step: 


- 52) = (El + E2 + Bsf - (Bi + B2 + Bg)^ 

= {E2 sin72 + E3 cos/^s)^ + (Essin/^s - Ei sin 71)^ + (Ei cos 71 + E2 cos 72)^ 

- (E2COS72 - Easin/^s)^ - (Eacos/^s - Ei cos 71)^ - (Ei sin 71 +E2sin72)^ 
= El + E| + E| — El — E| — E| + 2E1E2 (cos 7i cos 72 — sin 71 sin 72) 

+2E1E3 (cos 71 cos / 3 a — sin71 sin ( 3 ^) + 2E2Ea (cos 72 sin+ sin72 cos / 3 a) 

= 2E1E2 cos (71 + 72) + 2E1E3 cos (71 + Ps) + 2E2E3 sin (72 + P3) (B.l) 


(E • B) = (El + E2 + E3) • (Bi + B2 + B3) 

= (E2 sin 72 + Ea cos P3) (E2 cos 72 - E3 sin P3) 

+ (-El sin71 + Ea sin333) (-Ei cos71 + E3 cos/33) 

+ (El cos 7i + E2 cos 72) (-El sin 71 - E2 sin 72) 

= E1E2 (— COS71 sin 72 — cos 72 sin7i) 

+E1E3 (— sin71 cos Ps — sin/33 cos 71) 

+E2E3 (- sin 72 sin P3 + cos P3 cos 72) 

= -E1E2 sin (71 + 72) - E1E3 sin (71 + P3) + E2E3 cos (72 + /33) (B.2) 
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(^E^- E = [cos/^s cos (71 + 72) + sin72 cos (71 + /^a)] e*('^4-r-a;4t)^ 

+ E1E2E3 _ g-^^^ g-^ g^(k4•r-a;4t)^ 

+ Ei-E 2 - 5'3 _l_ _|_ _l_ gi(k 4 -r-i^ 4 t )2 (B_ 3 ) 

(^E^- B = J_ sin/Jg cos (71 + 72) + cos 72 cos (71 + P3)] e*(^ 4 -r-a; 4 t)^ 

+ E 1 -E 2 - 5'3 J^^g ^^g _|_ _ ^^g g^^^ _|_ ^ l ( k 4 ^ r ■ l ^4 t ) ^ 

+ Ei-E 2 - 5'3 j_ PQg gin (-^2 + Ps)] e*('^ 4 -r-a; 4 t) - 

(E • B) E = ^Qg gij^ + 72) - sin 72 sin (71 + /Jg)] e*('^4'r-a;4t)^ 

O 

+ Ei.E 2£^3 gin/ 3 g sin (71 + 72) - sin 71 cos (72 +/^s)] e*(^^'’'“‘^ 4 ‘)y 

O 

_ EiE2E^ + /^a) + cos 71 cos (72 + /^a)] z (B. 5 ) 


(E • B) B = [gin i^g sin (7^ + -)/2) - cos 72 sin (71 + / 3 g)] x 

O 

H—^ [-cos/?3 sin (71 + 72) - cos 71 cos (72 + Ps)] e*^ 

O 

+ EiE 2-£^3 [gin^2 sin (71 + P3) - sin 71 cos (72 + Ps)] (B.6) 


V- [(e^-b^ie 


EiE2E^ . 


ifc [(cos/^a + sin/Jg) cos (71 + 72) 

+ ^sin 72 - ^ cos 72^ cos (71 + P3) 


+ ( - sin 71 - - COS71 ) sin (72 + P3) 


gi(k4-r-(^4i) 


A 

dt 


(e^-bA e 

3 , 

= — ziv 

(e^-bA E 

Lv J \ 

2 

Lv J 1 


(B.8) 
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V X 


V X [(E-B)E] = 


V- [(E-B)B] 




+ 


4 

+ 2 /^3 cos (71 + 72) - - cos 7i sin (72 + /^s) 

, ,1 


gi(k4-r-a;4t) 


X 


+ 


E 1 E 2 -E 3 . 


^ sin /?3 cos (71 + 72) - ^ cos 72 cos (71 + f3s) 
+ sin 72 cos (71 + Ps) + sin 71 sin (72 + /^s)] 


■I* 

—ik [sin/Js cos (71 + 72) - cos 72 cos (71 + P^) 
+ coscos (71 +72) 

- COS71 sin (72 + /Ja)] 


(B. 9 ) 


-- —-ik [- cos 72 sin (71 + /^a) + cos 71 cos (72 + / 3 a) 

O 

sin/^a sin (71 + 72) - ^ sin 71 cos (72 + / 3 a) 
+ i sin 72 sin (71 +/ 3 a) 


gi(k4-r-aJ4i) 


X 


+ 


E 1 E 2 -E 3 . 


■ik 


-cos/ 3 a sin (71+72) , ^ 

(71 + P^) - cos 71 cos (72 + / 3 a)] e*(*^4'r-(^4t)^ 


+ 


+ cos 72 sin 

fp fp fp^ 

- -ik [cos / 3 a sin (71 + 72) + sin 72 sin (71 + / 3 a) 

-sin/ 3 asin (71 +72) 


- sin7i cos (72 +/ 3 a)] 


(B.IO) 


E1E2EI . 


ik [(sin / 3 a - cos / 3 a) sin (71 + 72) 

+ cos 72 - ^ sin 72^ sin (71 + / 3 a) 


+ ( - cos 71 + - sin 7i ) cos (72 + / 3 a) 


2 *(k 4 -r-‘-’ 4 t) (B.ll) 


^[(E.B)B] = -^ia;[(E.B)B] 


(B. 12 ) 
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(iVl) = 


X • {N2U) 


y • (iV234) 


_EiE2E^ . 


ik [(cos /^s + sin ( 5 ^) cos (71 + 72) 

+ ^sin72 - ^ cos 72^ cos (71 + ( 3 ^) 

+ sin 71 - ^ cos 71^ sin (72 + ( 3 ^) 

7 

+- (sin ( 3 ^ - cos ( 3 ‘i) sin (71 + 72) 

+ ^ (- cos 72 - ^ sin 72^ sin (71 + /^s) 
(- cos 71 + ^ sin71^ cos (72 + /^s) 


7(k4-r-a;4i) (B73) 




-- sin72 cos (71 + /^s) + sin71 sin (72 + ( 3 ^) 
-2 cos ( 3 ^ cos (71 + 72) + ^ cos 71 sin (72 + /Ss) 


7 

+ 4 


- cos 72 sin (71 + ^3)+ cos 71 cos (72 + Ps) 
-2sin/33 sin (71 + 72) 


sin 71 cos {j2 + P3) 


j(k4-r-u}4t) 


(B.14) 




-2 sin Ps cos (71 + 72) + ^ ^os 72 cos (71 + P3) 


- sin 72 cos (71 + A) + 2 sin 7 i sin (72 + Ps) 


+ - 


2 cos Ps sin (71 + 72) + - sin72 sin (71 + P3) 
+ cos 72 sin (71 + Ps) 
cos 71 cos {^2 + Ps) 


j(k4,-r-u}4t) 


(B.15) 
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z • {N 2 U) = i 


E1E2EI . 


ik { [- sin/33 cos (71 + 72) - cos ( 5 ^ cos (71 + 72) 


cos 72 cos (71 + /33) - ^ cos 71 sin (72 + ( 3 ^) 
7 

+ - [cos /33 sin (71 + 72) - sin/33 sin (71 + 72) 


--sin 72 sin (71 +/33) 
+ isin7i cos (72 + Ps) 


7(k4-r-aj4t) 


(B. 16 ) 


V(/V 1 ) =ik{Nl) 


(x + y 



(B. 17 ) 


o Q 

-(iv^ = --iupmuj 


(B.18) 


By plugging these expressions into Eq. ( 3 . 38 ), we finally arrive at Eq. ( 5 . 19 ). 
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Appendix C 


Geometric factors in 
hep-ph/0510076 


Below follow the geometric factors refered to in the paper Using high-power lasers for de¬ 
tection of elastic photon-photon scattering (hep-ph/ 0510076 ). Due to somewhat different 
conventions, Gs^ here differs by a constant factor from the expression presented above in 
the MSc thesis. 

For wave vectors confined to the x-y-plane, kj=cos(/>jX + sinc/j-y, and with 7^ defined 
as the angle between the z-axis and the polarization vector of the Sj-field, in such a way 
it is positive in clockwise direction when looking in the kj-direction, the two dimensional 
geometric factor becomes 


G2d - +2 


7 

- cos 73 sm 


+ 


(71 + 72) - sin73 cos (71 -b 72) 
[sin ((/)3 — (f>4) cos (1)4 — (sin (j)^ — sin (j)4)] sin^ 
in (71 -I- 73) - sin72 cos (71 + 73) 


7 

- cos 72 sm 


X sin (02 - 


+ 


(02 — 04 ) COS 04 — (sin 02 — sin 04 )] sin^ 
;in (72 -I- 73) - sin71 cos (72 + 73) 


- cos 71 sin 


+ - 


+ 


+ 


[sin (01 — 04) cos 04 — (sin 0i — sin 04)] si 
7 

- cos 73 sin (71 + 72) - sin 73 cos (71 -|- 72) 
[sin (03 — 04) sin 04 -|- (cos 03 — cos 04)] sin^ 
^ cos 72 sin (71 + 73) - sin 72 cos (71 + 73) 
[sin (02 — 04) sin 04 -)- (cos 02 — cos 04)] sin^ 
^ cos 7i sin (72 + 73) - sin 71 cos (72 + 73) 


01 — 02 


01 - 03 


02 — 03 


01 — 02 


01 - 03 


X 
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X [sin ((/)! — ^4) sin (j)^ + (cos (pi — cos (p^)] sin^ 

I ^ sin 73 sin (71 + 72) + cos 73 cos (71 + 72) 

■ 2 ^ 03 - </>4 ^ ^.^2 ^ 01 ~ 02 ^ 


^ 02 — 03 ' 


X sin^ 


+ 


7 

- sin 72 sin (71 + 73) + cos 72 cos (71 + 73) 
^ 02 - 04 \ ^.^2 f 01 ~ 03 \ 


X sin^ 


+ 


■7 

- sin 71 sin (72 + 73) + cos 71 cos (72 + 73) 
. 2 02 “ 03 ^ 'I 

Sin - 


X sin 


^ 01 — 04 " 


• Z. 


(C.l) 


With wave vectors ki = /cx, k 2 = ky, k 3 = |z, k 4 = k'k + ky — ^kz, and 03 defined 
as the angle between the x-axis and the polarization vector of the £' 3 -field, in such way 
that it is positive in clockwise direction when looking in the k 3 -direction, we get the three 
dimensional geometric factor 

r - ^ 

- - 7 : 


9 


^ sin 03 - cos 03 ^ cos (71 + 72) 

1 


i cos 72^ cos (71 +03) 


+ - sin 72 - T cos 


7 
+ 4 


2 

'9 


+ Q sin 71 + ^ cos 71^ sin (72 + 03) 

7 r/ 1 \ 

- sin 03 - 2 CO® j sin (71 + 72) 

cos72 - ^ sin72^ sin (71 + 03) 

Q cos 7i - ^ sin71^ cos (72 + 03) 

1 \ 


+ 

+ 


> X 


1 - sin 03^ cos (71 +72) 


i cos 03 

+ ("i cos 72 - ^ sin72cos (71 + 03 ) 


+ 


7 

+ 4 


+ 


9 


\ a. / 

COS 7 i - ^ sin71^ sin (72 + 03 ) 
^cos 03 + ^ sin 03^ sin (71 + 72) 

+ Q sin 72 + ^ cos 72^ sin (71 + 03) 
+ Q sin 7i - ^ cos 71^ cos (72 + 03) 
^ (cos 03 + sin 03) cos (71 + 72) 

+ i (sin 72 + cos 72) cos (71 + 03) 

O 
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+- (cos7i - sin7i) sin (72 + /^s) 


7 

+ 4 


- (sin /Js - cos /^s) sin (71 + 72) 

+ ^ (sin72 - cos 72) sin (71 + (3^) 

O 

+ ^ (- cos 71 - sin 71) cos (72 + /^s) 


(C.2) 


For this configuration the number of scattered photons per shot is estimated to 


= 1.31r/'Gi, 


1 ;Um 


L 


1 //m 


PlP2P^ \ 

1PW3 ) 


(C.3) 


where Pj is the power of the incoming pulses, L the pulse length, A 4 the generated wave¬ 
length, and with rf given by 


2 _ 72 - r lg{r,e,4>) 

Jo Jo Ig^resix^ 
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(•27r /-TT sm 


sin 0 d6 d(j) 

/C 4 I (1 — cos(/>sin 0 ) 


sm 


/C 4 I (sin (j) sin 6) 


sm' 


0 Jo (1 — cos0sin0)^ 

2 '/C 4 I (cos 0)1 


(sin 0 sin 0)2 


(cos0)" 


sin0 dO dcj), 


(C.4) 


where Ig is the generated intensity, Ig^res its maximum value in the resonant direction, 
and b is the pulse width. 

An optimal choice of polarization angles for the three dimensional configuration is 
given by 71 = 0, 72 = = 7r/2, yielding Gl^ = 0.77. 
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